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VAN  VRAAG  NAAR  ANTWOORD 


Met  hehulp  van  gevechtssimulaties  kunnen  tal  van  viagen  bcantwoord  wurden.  Vooibeelden  van 
zulke  vragen  zijn; 

Wat  is  de  invloed  van  een  bepaalde  vcrbetering  aan  een  wapensysteem? 

Welke  wijze  van  optreden  is  bcter? 

Welke  mix  van  wapensystemen  is  het  bcste? 

In  dit  hoofdstuk  wordt  beschreven  wat  er  gedaan  meet  worden  om  met  behulp  van  een 
stochastiscb  simulatiemofiel  antwoord  tc  Vrijgcn  op  bovenstaande  vragen,  als  al  in  e-e.:  ecrder 
stadium  is  gekozen  voor  een  stochasdsch  simulademodel  voor  het  haantwoorden  van  de  vraag. 
Aan  het  eind  van  het  hoofdstuk  is  de  beschreven  methodiek  wecrgegeven  in  twee  schema's.  Het 
eerste  schema  (figuur  1)  geeft  weer  hoe  tot  een  keuze  van  een  stadstische  methode  kan  worden 
gekomen.  Het  twecde  schema  (figuur  2)  geeft  weer  hoe  onderzocht  kan  worden  of  het 
waarschijnlijk  is  dat  een  set  waamemingen  uit  een  noimaJe  verdeling  komt. 

Om  te  beginnen  moet  de  vraag  duidelijk  zijn.  Een  aantal  dingen  moet  bij  de  vraagstelling  bepaald 
worden: 

Op  basis  van  welke  maat  van  effecdviteit  (MOE)  bepaalde  uitkomsten  met  elkaar  worden 
vergeleken. 

Op  welk  moment  de  MOE  uitgerekend  moet  worden. 

Dit  kan  het  einOe  van  de  simulade  zijn  of  bepaald  worden  aan  de  band  van  een  aantal 
criteria. 

In  cen  gevecMssimulatie  kan  bijvoorbeeld  de  MOE  uitgerekend  worden  als  de  ene  partij  tot 
40%  gesleten  is  of  als  de  andere  paitij  tot  30%  gesleten  is. 

Vaak  kan  er  geen  antw  .,/rd  gekrcgen  worden  met  belulp  van  de  bestaande  scenario's.  Met  een 
scenario  wordt  bedoeld  de  hele  set  van  invoergegevens  die  nodig  is  om  een  simulatie  te  draaien. 

Er  zuUeu  nieuwe  scenario's  gecreiierd  moeten  worden  of  oude  aangepast.  Voor  het  beantwoorden 
van  de  meeste  vragen  zijn  mcerdere  scenario's  nodig.  Door  het  vergelijken  van  de  uitkomsten  van 
de  scenario's  kunnen  aiuwoorden  verkregen  worden. 
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"or  answering  C(uestions  about  the  deployment  of  weapotj^ystems,  the  stochastic  combat  simulation  model 
FSM  (Force  Structure  Model)  is  one  of  the  models  developed.  This  report  describes  the  way  to  answer  a 
question  by  means  of  the  stochastic  simulation  model  and  statistical  techniques. 

The  outcomes  of  combat  simulations  have  to  be  valued  before  they  can  be  compared  to  each  other.  This  can 
be  done  with  the  help  of  one  or  more  measures  of  effectiveness  (MOE’s).  When  using  stochastic  model,  the 
model  has  to  be  runned  several  times  with  the  same  input.  In  this  way  a  set  of  observations  is  generated  and 
something  can  be  said  about  the  distribution  of  possible  outcomes. 

In  this  'cpnrt  a  number  of  statistical  tcots  are  described  with  which  two  sets  of  observations  can  be  compared 
to  car., it  other. 

Some  computer  programs  are  written  for  making  the  analysis  of  simulation  results  easier.  These  programs 
are  also  described. 

An  impulse  is  given  to  further  research.  — — 
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1  INLEIDING 

Binnen  de  researchgroep  wapensysteem-  cn  beleidsstudies  LAND  van  hct  Fysisch  cn 
Elektronisch  Laboratorium  wordt,  voor  het  uitvocren  van  studies,  gcbruik  gcmaakt  van 
stochastiscbe  gevechtssimulatiemodellen.  Het  Force  Stnicture  Model  (FSM)  is  66n  van  die 
modellen. 

In  het  kader  van  de  opdracht  FSM  (A89KL619)  is  onderzoek  gedaan  naar  de  metbodiek  om  met 
behulp  van  dergelijke  modellen  antwoorden  te  verkrijgen  op  vragen  die  voortvloeien  uit  de 
studies.  Deze  metbodiek  wordt  in  hoofdstuk  2  beschreven. 

De  metbodiek  is  gebasecrd  op  de  kwantitatieve  waardering  van  de  simolatieresultaten.  Dit  houdt 
in  dat  eUce  simulatie  uitkomst  in  ecn  getal  wordt  uitgedrukt.  Zo'n  getal  wordt  MOE  (Measure  Of 
Effectiveness)  genoemd.  In  hoofdstuk  3  worden  enkele  MOE's  voor  gevechtssimulatief 
beschreven. 

Met  behulp  van  MOE's  is  het  mogelijk  om  simulatieresultaten  grafisch  wecr  te  geven  en  met 
elkaar  te  vergelijken.  Hierbij  wordt  gebruik  gemaakt  van  statistische  methoden.  Deze  metboden 
worden  in  hoofdstuk  4  beschreven. 

Om  de  analyse  van  simulatieresultaten  te  vergemakkelijken  zijn  drie  computerprogramma's 
geschreven.  In  hoofdstuk  5  worden  deze  drie  programma's  beschreven.  Ecn  aanzet  tot 
vervolgonderzoek  wordt  gedaan  in  hoofdstuk  6. 

De  hoofdstukken  zijn  afzonderlijk  te  lezen.  Aihankelijk  van  de  inieresse  van  de  lezer  kan  bet 
lezen  van  een  enkel  bonfdstuk  voldormde  zijn. 

Via  de  index  in  bijlage  A  kan  gevonden  worden  waar  in  bet  rapport  meer  over  een  bepaald 
onderwerp  staat  geschreven. 
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Reladeve  gevechtskracht  blauw 

De  reladeve  gevechtskracht  blauw  woidt  gedefinieerd  als  de  huidige  gevechtskiacht  van  blauw 
gedeeld  door  de  initiele  gevechtskracht  van  blauw.  Deze  MOE  geeft  alleen  de  eigen  sterkte  weer. 
In  het  vervolg  wordt  deze  MOE  aangeduid  met  FORCE_BLUE.  FORCE_BLUE  kan  als  voIgt 
worden  uitgedrukt: 


B-b 

FORCE  BLUE  =  — 

B 

Het  bereik  van  deze  MOE  is  [0,1]. 

Relatieve  verliezen  rood 

Reladeve  verliezen  rood  wordt  gedefinieerd  als  de  verliezen  van  rood  gedeeld  door  de  inidble 
gevechtskracht  van  rood.  Deze  MOE  geeft  het  verlies  van  de  vijand  weer,  en  kan  als  volgt  worden 
uitgedrukt: 


LOSS_RED=  - 
R 

Ook  deze  MOE  ligt  tussen  0  en  1. 

Reladeve  gevechtskracht  verfaouding 

De  reladeve  gevechtskracht  verbouding  wordt  gedeEnieerd  als  de  verbouding  van  de 
gevechtskracht  van  blauw  tot  de  gevechtskracht  van  rood  gedeeld  dow  de  inidSle  gevechtskracht 
verbouding.  Deze  MOE  geeft  de  reladeve  stedrte  van  blauw  ten  opzicbte  van  rood  weer.  In 
formulevoim  ziet  deze  MOE  er  als  volgt  uit: 


FORCE_RATIO  = 


B-b/R-r 

"b^ 


Deze  MOE  ligt  tussen  0  en  oneindig. 
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Met  elk  scenario  moet  een  aantal  runs  gediaaid  worden.  Het  aantal  luns  bemvloedt  de 
betrouwbaarheid  van  de  in  het  stadstisch  deel  gebniikte  schattingen.  Het  diaaien  van  runs  kost 
echter  veel  (computer)tijd  en  veel  geheugencapaciteit  Over  het  algemeen  zuUen  50  runs  een 
goede  indicatie  geven  over  de  verdeling  van  de  uitkomsten. 

Bereken  nu  de  MOE  van  elke  run. 

Maak  van  de  MOE's  voor  elk  scenario  een  frequentiehistogram  en/of  een  grariek  van  de 
empirische  verdelingsfunctie.  Dit  geeft  een  indicatie  of  er  lets  vreemds  aan  de  hand  is  (rare 
uitschietrrs  en  dergelijke  kunnen  zo  opgespoord  worden). 

Als  er  met  verschillende  criteria  is  gewerkt  voor  het  bepalen  van  het  moment  waarop  de  MOE 
uitgerekend  moest  worden  is  het  zinvol  om  fiequentiehistogrammen  te  maken  die  bij  elke  run 
aangeven  om  welke  reden  de  MOE  van  die  run  op  dat  moment  is  uitgerekend. 

Nu  kunnen  de  scenario's  paarsgewijs  vergeleken  worden.  Bij  elke  vergelijking  hoort  een 
nulhypothese,  een  altematieve  hypotbese  (Meestal:  nulhypothese  is  niet  waar)  en  een 
betrouwbaarheidsniveau  1  •  a. 

a  geeft  aan  wat  de  toelaatbare  kans  is  dat  de  altematieve  hypothese  aanvaard  wordt  terwijl  de 
nulhypothese  waar  is.  0.05  en  0.10  zijn  gehraikelijke  waarden  voor  a. 

Er  kunnen  twee  verschillende  groepen  nulhypotheses  onderscheiden  worden.  Ten  eerste  de  groep 
waarbij  de  hele  verdeling  van  de  uitkomsten  van  belang  is.  Dan  luidt  de  nulhypothese 
bijvoorbeeld:  de  verdeling  van  de  MOE's  van  het  ene  scenario  is  gelijk  aan  de  verdeling  van  de 
MOE's  van  het  andere  scenario.  Ten  tweede  de  groep  waarbij  alleen  de  verwachting  van  belang 
is.  Dan  luidt  de  nulhypothese  bijvoorbeeld:  de  verwachting  van  de  MOE  in  bet  ene  scenario  is 
gelijk  aan  het  verwachting  van  de  MOE  in  het  andere  scenario. 

Het  vergelijken  van  de  gehele  verdelingen  van  de  MOE's  uit  twee  scenario's  begint  met  het  maken 
van  tea  frequentiehistogram  en  een  gra&ek  van  de  empirische  verdelingsfuncties.  Op  het  oog  kan 
dan  al  ieu  gezegd  worden  over  de  overeenkomst  tussen  de  verdelingen  van  de  MOE's.  Dit  is 
echter  zeer  subjectief.  Daarom  moet  er  daama  de  toets  van  KolmogoroV'Smimov  gedaan  worden. 
In  het  programma  DEC_M0E_V1S  is  een  gewijzigde  KolmogoroV'Smimov  toets  ingebouwd. 
Deze  toets  weet  raad  met  nulhypothetes  als:  de  verdeling  van  de  MOEs  van  het  ene  scenario 
verschilt  minder  dan  een  door  de  vraagsteller  vastgestelde  waarde  van  de  verdeling  van  de  MOE's 
uit  bet  andere  scenario. 
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Om  alleen  de  verwachtingen  met  elkaar  te  vergelijken  moet  east  een  stadstisch  overzicht 
gemaakt  worden  van  de  MOE's  uit  beide  scenario's.  Hieruit  kunnen  de  gemiddelde  worden 
afgelezen.  Dit  zijn  schatters  voor  de  verwachting.  Het  vaschil  mssen  de  twee  gemiddelden  geeft 
een  indicade  voor  het  al  dan  niet  aanvaarden  van  de  aanname  dat  a  verschil  is  tussen  beide 
verwachdngen. 

Daama  moet  bekeken  worden  of  de  waamemingen  van  de  MOE's  afkomsdg  kunnen  zijn  uit  een 
normale  verdeling.  Dit  is  bepalend  voor  het  latae  gebruik  van  toetsen  die  gebaseerd  zijn  op  de 
normale  vadeling  of  voor  het  gebruik  van  vadelingsvrije  toetsen. 

Het  testen  of  de  waamemingen  afkomsdg  kunnen  zijn  uit  een  nonnale  vadeling  begint  aan  de 
hand  van  graflsche  weagaven  van  de  resultaten  en  aan  de  hand  van  de  stadsdeken  scheefheid 
(Engels:  skewness)  en  kuitosis.  De  scheefheid  en  kurtosis  zijn  gelijk  aan  0  voa  een  normale 
verdeling.  De  beide  stadsdeken  mogen  niet  al  te  veel  van  0  afwijken  om  aan  te  nemen  dat  de 
waamemingen  uit  een  nomale  verdeling  komen.  Voor  SO  waamemingen  zou  als  richtlijn  een 
maximale  afwijidng  van  1.0  genomen  kunnen  worden. 

Als  nu  nog  steeds  getwijfeld  wordt  ova  het  al  dan  niet  afkorisdg  zijn  van  de  waamemingen  uit 
een  normale  vadeling  dan  moet  een  Goodness  of  fit  toets  uitgevoad  worden.  Als  a  veel 
waamemingen  zijn  dan  kan  het  beste  de  Chi-kwadraat  toets  gebruikt  worden  (richtlijn:  mea  dan 
25  waamemingen).  Zijn  a  minda  waamemingen  dan  is  het  beta  om  dn  Kolmogorov-Smimov 
toets  te  gebruiken.  Het  gebruik  van  deze  toets  is  echter  af  te  raden  als  de  waamemingf'  jit  ees 
disaete  vadeling  lijken  te  komen.  dat  geval  moet  de  aanname  dat  de  waamemingen  uit  een 
nonnale  vadeling  komen  verwotpen  waden. 

In  Hguur  2  staat  de  methodiek  schematisch  weergegeven. 

Als  in  een  eerda  stadium  de  aannames  zijn  aanvaard  dat  de  waamemingen  van  MOE's  uit 
normale  vadelingen  komen,  kan  een  Student's  t-toets  gebruikt  worden  of  een  variant  hierop,  de 
Satterthwate  toets.  Deze  laatste  mag  gebruikt  worden  zonda  de,  op  tbeorediche  grooden 
gebaseade,  aanname  van  gelijke  varianties. 

Is  de  aanname  van  nattule  verdelingen  niei  gerecbtvaardigd  dan  moet  een  vadelingsvrije 
metbode  gebruikt  worden.  Dit  kan  zijn  de  toeu  van  Wilcoxon  of  de  bootstrap  inohode. 

Met  bchulp  van  de  toets  van  Wilcoxon  kan  getoetst  worden  of  de  vowacbtin^o  van  elkaar 
vetschUlen  en  met  behulp  van  de  boouhap  metbode  kan  een  betrouwbaarfaeidsinterval  geaeSerd 
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Worden  vou  c';i  verschil  tussen  de  verwachtingen.  Dit  interval  geeft  aan  dat  met  de  opgegeven 
leans  de  werkelijke  waarde  van  het  verschil  tussen  de  vowachtingen  in  het  interval  ligt. 

Soms  lean  het  nodig  zijn  dat  van  6in  of  meerdere  scenario's  een  aantal  extra  runs  gedraaid  moeten 
Worden  om  een  antwoord  met  meer  precisie  te  veikiijgen. 

Nu  moet  met  behulp  van  de  veikregen  stadstische  gegevens  antwoord  worden  verkregen  op  de 
vraag.  Het  is  moeilijk  om  in  het  algemeen  te  beschrijven  hoe  dit  moet.  In  het  geval  van  een 
eenduidige  vraag  (bijvoorbeeld:  duurt  het  in  dit  ene  scenario  net  zolang  tot  de  vijand  is  gesleten 
tot  40%  dan  in  dat  andere  scenario)  dan  is  een  antwoord  vrij  makkelijk  te  given.  Vaak  is  de 
vraagstelling  complexer  en  roept  een  stadstische  vergelijking  andere  vragen  op.  Ook  is  het  vaak 
zo  dat  er  niet  vergeleken  wordt  op  basis  van  ^n  MOB  maar  op  basis  van  meerdere.  In  dat  geval 
moeten  technieken  als  multi-criteria  analyse  uitkomst  brengen.  Stadstische  technieken  kunnen 
dan  gebraikt  worden  om  interval-schattingen  te  maken  voor  verwachdngen  of  voor  verschillen 
tussen  verwachtingen  van  twee  scenario's. 
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3  MATEN  VAN  EFFECIIVITEIT 

Dit  hoofdstuk  gecft  ecn  beschrijving  van  verschiUende  maten  van  effectiviteit  (MOE's)  voor 
gevechtssimulaties.  Eerst  wordt  ecn  overzicht  van  ecn  aantal  MOE’s  gegcven.  Vervolgens  worden 
dcze  MOE's  nader  beschrevcn,  en  wordt  ei  langegeven  in  welic  soort  gevecht  ze  gebruikt  kunnen 
worden.  De  geschiktheid  van  een  MOE  hangt  namelijk  af  van  bet  soort  gevecht  dat  gevo.xd 
wordt. 

3.1  Overzicht 

Enkele  veelgebruikte  MOE's  zijn; 

1.  Relatieve  gevechtskracht  blauw 

2.  Relatieve  verliezen  rood 

3.  Relatieve  gevechtskracht  verhouding 

4.  Relatieve  verlies  verhouding 

5.  Win  factor 

6.  Succjs  boolean 

7.  TerreinwLnst 

8.  Tijdsduur 

Door  ODS  zclf  werd  ontwikkcld: 

9.  Hindemis  verdedigingskracht 

Deze  MOE's  worden  in  de  volgende  paragraaf  uitgewerict. 

3.2  Nadere  beschrijving  van  de  MOE's 

Ont  te  komen  tot  een  duidelijke  beschrijving  van  de  in  3.1  genoemde  MOE's,  wnden  de  volgende 
definities  gebruikt: 

B  ;  Inidele  gevechtskracht  blauw 

R  ;  Inidele  gevechtskracht  rood 

b  :  Verlies  blauw 

r  :  Verlies  rood 

waarbij  de  gevechtskracht  en  de  verliezen  worden  uitgediukt  in  aantalkn  voertulgen.  Met  behulp 
van  deze  definides  worden  nu  de  verschiUende  MOEs  beschrevcn. 
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Relatieve  eevechtskracht  blauw 

De  relatieve  gevechtskracht  blauw  wordt  gedefinieerd  als  de  huidige  gevechtskracht  van  blauw 
gedeeld  door  de  initiele  gevechtskhicht  van  blauw.  Deze  MOE  geeft  alleen  de  eigen  sterkte  weer. 
In  het  vervolg  wordt  deze  MOE  aangeduid  met  FORCE_BLUE.  FORCE_BLUE  kan  als  volgt 
worden  uitgedrukt: 


B-b 

FORCE_BLU2=  — 

B 

Het  bereik  van  deze  MOE  is  [0,1]. 

Relatieve  verliezen  rood 

Relatieve  verliezen  rood  wordt  gedefinieerd  als  de  verliezen  van  rood  gedeeld  door  de  initiiile 
gevechtskracht  van  rood.  Deze  MOE  geefl  het  verlies  van  de  vijand  weer,  en  kan  als  volgt  worden 
uitgedrukt; 


LOSS_RED=  - 
R 

Ook  deze  MOE  iigt  tussen  0  en  1. 

R£larissa.gty.«hi^Mbt  .vsriwudag 

De  relatieve  gevechtskracht  verhouding  wordt  gedefinieerd  als  de  veibouding  van  de 
gevechtskracht  van  blauw  tot  de  gevechtskracht  van  rood  gedeeld  door  de  initi£le  gevechtskracht 
verhouding.  Deze  MOE  geeft  de  relatieve  stedUe  van  blauw  ten  opzichte  van  rood  weer.  In 
formulevorm  ziet  deze  MOE  er  als  volgt  uit: 


FORCE_RAT10  = 


B-b/R-r 

b/r' 


Deze  MOE  ligt  tussen  0  en  oneindig. 
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Relatieve  verlies  verhoudiag 

De  felatievc  verlies  verhouding  wordt  gedefinieerd  als  de  verhouding  van  de  verliezen  van  rood 
tot  de  verliezen  van  blauw  gedeeld  door  de  inidele  gevechtskracht  verhouding.  Deze  MOE  geeft 
de  relatieve  verliezen  van  rood  ten  opzichte  van  blauw  weer,  en  kan  als  volgt  worden  uitgedrukt: 


r  h 

LOSS_RATIO=  - 

R/B 

Deze  MOE  ligt  tussen  0  en  oneindig. 

Win  factor 

De  win  factor  is  een  combinatie  van  de  relatieve  gevechtskracht  verhouding  en  de  relatieve 
verlies  verhouding.  Deze  kan  als  volgt  geschreven  worden; 


WIN_FACTOR  = 


FORCE.RATIO  +  LOSS.RATIO 
2 


Ook  de.  i  MOE  ligt  tussen  0  en  oneindig. 

Opdracht  geslaayd 

Deze  MOE  geeft  weer  of  een  opdracht  geslaagd  is  of  niet.  Wanneer  de  (^)dracht  geslaagd  is,  krijgt 
deze  MOE  de  waarde  1;  wanneer  de  opdracht  niet  geslaagd  is,  de  waarde  0. 

SUCCES=  f  ^  als  opdracht  geslaagd 

I  0  als  opdracht  niet  geslaagd 

Deze  MOE  kan  dus  alleen  de  waarde  0  of  1  aannenoen. 

iHKiflwinsl 

Deze  MOE  geeft  weer  hoeveel  reireinwinst  er  door  de  aanvallende  party  geboekt  is.  Ervan 
uitgaande  dat  rood  de  aanvallende  paitij  is,  kan  er 
geschreven  worden: 
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TERRAIN  =  min  (x-coordinaten  rood) 

Bij  de  bepaling  van  deze  MOE  wordt  aangenomen  dat  rood  uit  het  oosten  aanvalt.  Deze  MOE 
wordt  groter  naarmate  rood  meer  terreinwinst  boekt. 

Is  blauw  de  aanvallende  partij,  dan  kan  de  terreinwinst  als  volgt  uitgedrukt  worden 


1 

\ 

i 


TERRAIN  =  max  (x-cooidinaten  blauw) 


waarbij  ervan  uitgegaan  wordt  dat  blauw  uit  bet  westen  aanvalt. 


Tijdsduur 

De  tijdsduur  geeft  weer  hoe  lang  blauw  zijn  verdediging  beeft  gebandbaafd.  Wanneer  rood  zich 
tenigtrekt  beeft  deze  MOE  een  vaste  waarde  T,  wanneer  blauw  zich  terugtrekt  geeft  deze  MOE 
het  tijdstip  t  waarop  blauw  zich  terugtrekt: 


DURATION  = 


T 

t 


als  rood  zich  terugtrekt 
als  blauw  zich  terugtrekt 


waarin  T  constant  en  T  >  t. 


Hindemis  verdedigingskracht 

De  hindemis  verdedigingskracht  geeft  een  waardering  aan  een  opstelling  van  eenbeden,  die  een 
hindemis  verdedigen.  Hiertoe  wordt  berekend  wat  de  gemiddelde  kana  is,  dat  eM  vijandelijke 
tank  die  de  hindemis  betreedt,  wordt  uitgeschakeld.  Hierbij  wordt  in  beschouwing  genomen 
welke  verdedigende  eenbeden  zicht  hebben  op  (delen  van)  de  hindemis. 

De  berekening  van  deze  kans  geschiedt  door  de  hindemis  op  te  delen  in  N  stukjes  van  S  meter. 
Vcrvolgens  wordt  voor  elk  stukje  bj  bepaald  welke  blauwe  eenbeden  liue^f-si^t  met  dit  stukje 
hindemis  hebben.  Voor  elke  eenheid  j  met  line-of-sight  wordt  bepaald  wid  de  kill>kant  Pg  is 
tegencver  een  tank  die  zich  op  hi  bevindt 
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Met  behulp  van  de  fomule 


Pi=i-  na-Py) 

j 

kan  nu  de  kill-kans  tegenover  een  tank  op  h|,  worden  beiekend.  De  gemiddelde  Idll-kans  wordt  nu 
berekend  door  het  gemiddelde  over  alle  stukjes  hj  te  nemen; 


1 

OBST  DEFENCE  =-  IP 

N  1 

Deze  MOE  ligt  tussen  0  en  1. 

3.3  Bruikbaarheid  van  de  MOE's 

De  bruikbaarheid  van  de  MOE's  hangt  af  van  het  sooit  gevecht  dat  gevoerd  wordt.  De  volgende 
gevechtsvormen  worden  onderscheiden: 

1.  Verdedigend  gevecht: 

Gevechtsvorm  die  ten  doel  heeft  de  aanvaller  zo  mogelijk  de  toegang  tot,  maar  tenminste  de 
doorgang  door  het  ter  verdediging  toegewezen  vak  -  onder  behoud  van  het  tactiach  essentieel 
gcbied  -  tc  ontzeggen. 

2.  Vertragend  gevecht: 

Gevechtsvorm  met  als  doel  de  aanvaller  over  een  opgedragen  diepte  van  het  toegewezen  vak 
op  te  houden  en  hem  verliezen  toe  te  brengen,  zodat  het  vervolgens  voeren  van  een 
verdedigend  of  aanvallend  gevecht  wordt  begunstigd. 

3.  Aanvallend  gevecht: 

Gevechtsvorm  waarbij  de  beschikbare  tniddelen  op  een  door  de  (hogere)  commandant  te 
bepalen  tijd  en  plaats  worden  ingezet  met  alt  doel  bet  vetmeesteren  van  een  tactisch 
essentieel  gebied,  waarbij  de  vijand  de  mogelijkheid  en  de  wil  wotdt  ontnomoa  verder 
weerstand  te  bieden. 
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Binnen  het  verdedigend  gevecht  bestaat  de  mogelijkheid  om  een  tegenaanval  of  een  tegenstoot  te 
plaatsen: 

Tegenaanval: 

Een  in  het  operadeplan  al  dan  niet  voorziene  aanval  met  het  doel  tactisch  belangrijk  gebied 
te  hememen  of  verloren  gaan  van  tactisch  essentieel  gebied  te  voorkomen. 

Tegenstoot: 

Een  in  het  operatieplan  al  dan  niet  voorziene  aanval  met  het  doel  de  vijand  zoveel  mogelijk 
verliezen  toe  te  brengen. 

De  eerste  vijf  MOE's  uit  het  overzicht  in  3.1  kunnen  bij  elk  sooit  gevecht  worden  gebniikt  De 
eigen  sterkte  en  de  sterkte  van  de  vijand  zijn  namelijk  aldjd  belangrijk,  ongeacht  het  soort  gevecht 
dat  gevoerd  wordt.  De  zesde  MOE,  succes  boolean,  lean  eveneens  in  elk  soort  gevecht  worden 
gebniikt,  aangezien  het  uitvoeren  van  een  opdracht  in  iedere  gevechtsvorm  voorkomt. 

Naast  deze  zes  algemene  MOE's  zijn  er  ook  MOE's  die  gebonden  zijn  aan  het  soort  gevecht  dat 
gevoerd  wordt.  Wordt  er  een  verdedigend  gevecht  gevoerd,  dan  is  de  terreinwinst  van  de  vijand 
belangrijk.  Hier  kan  dus  de  zevende  MOE,  terreinwinst,  worden  gebniikt.  Dient  er  in  het 
verdedigend  gevecht  een  hindemis  te  worden  verdedigd,  dan  is  de  hindemis  verdedigingskracht 
het  meest  geschikt  als  MOE. 

In  het  vertragend  gevecht  dient  de  vijand  zo  lang  mogelijk  opgehouden  te  worden.  De  tijdsduur  is 
hier  dus  een  geschikte  MOE.  Verder  dienen  er  zoveel  mogelijk  verliezen  aan  de  vijand  te  worden 
toegebracht.  Hierbij  zijn  de  reladeve  verliezen  rood  en  de  reladeve  verlies  verfaouding  gesclukte 
MOE's. 

Het  aanvallend  gevecht  heeft  als  doel  het  bezetten  van  vijandelijk  gebied.  De  terreinwint:  is  hier 
dus  een  geschikte  MOE.  Deze  MOE  kan  ook  worden  gebroikt  bij  de  tegenaanval. 

De  tegenstoot  heeft  als  doel  de  vijand  zoveel  mogelijk  verliezen  toe  te  brengen,  teneinde  de 
gevechtskrachtverhouding  in  bet  eigen  votnrdeel  1*  wijzigen.  MOE's  die  hierbij  gebniikt  kunnen 
worden  zijn  de  reladeve  verliezen  rood,  de  relatieve  verlies  vertiouding,  de  reladeve 
gevecbtskracht  verhouding  en  de  win  factor. 

In  sommige  situades  kan  bet  zinvol  zijn  om  een  combinade  v.in  MOE  's  te  gebiuiken. 
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4  GEBRUKTE  STAHSHSCHE  METHODEN 

In  dit  hoofdstuk  worden  enkele  methoden  bescbieven,  waairoee  twee  scenario's  op  basis  van  een 
gekozen  MOE  met  elkaar  vergeleken  kunnen  worden.  Hieitoe  woiden  voor  beide  scenario's  nit  de 
gedraaide  simulatieruns  de  waaiden  van  de  gekozen  MOE  bepaald.  Zijn  de  veikiegen  datasets 
normaal  veideeld,  dan  kunnen  verdelingsafhankelijke  methoden  worden  gebruikt.  Zijn  de 
verkregen  datasets  niet  normaal  verdeeld,  dan  dient  er  voor  een  verdelingsvrije  methode  te 
worden  gekozen. 

1 

j 

j  Cm  te  weten  te  komen  of  de  waamemingen  ncnnaal  verdeeld  zijn,  is  bet  zinvol  om  de 

I  waamemingen  grarisch  weer  te  geven.  Dit  kan  door  middel  van  frequentiehistogrammen  en 

1  empirische  verdelingsfuncties.  Vervolgens  kan  er  met  behulp  van  zogenaamde  "Goodness  of  rit" 

toetsen  worden  bepaald  of  het  waarschijnlijk  is  dat  de  waamemingen  uit  een  normale  verdeling 

} 

afkomstig  zijn. 

4.1  Frequendehistogrammen 

Voor  bet  maken  van  een  firequendehistogram  worden  de  waamemingen  in  een  aantal  klassen 
opgesplitst.  Het  aantal  waamemingen  binnen  een  klasse  noemt  men  frequende.  De  verdeling  van 
dc  fiequendes  over  de  verschillende  klassen  wordt  frequendeverdeling  genoemd.  Worden  de 
frequendes  gedeeld  door  het  totaal  aantal  waamemingen,  dan  spreekt  men  van  een  reladeve 
drequendeverdeling. 

Door  de  frequendeverdeling  uit  te  zetten  in  een  histogram  krijgt  men  inzicbt  in  de  werkelijke 
I  verdeling  van  de  waamemingen. 

j  Het  bovenstaande  wordt  ge'illustreerd  aan  de  band  van  bet  volgende  voorbeeld: 

Met  een  bepaald  scenario  werden  50  simuladeruns  gedraaid.  Uit  de  resultaten  hietvan  werd  voor 
elke  simuladenm  de  reladeve  gevechtskracbt  veriiouding  berekend.  De  SO  waamemingen  van 
deze  MOE  werden  in  6  klassen  verdeeld  waaibij  de  volgende  reladeve  frequendes  werden 
waargenomen: 
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FORCE_RATIO 

Frequentie  (relatief) 

0.5  <=1.0: 

11  (22.00%) 

1.0  <=  1.5: 

23  (  46.00%) 

1.5  <=  2.0: 

11  22.00% 

2.0  <=  2.5 : 

3  •  (  6.00% 

2.5  <=  3.0: 

0  0.00%) 

3.0  <=  3.5 : 

2  (  4.00%) 

In  de  volgende  figuur  zijo  de  relatieve  frequenties  uitgezet  in  een  histogram,  waarbij  tevens  de 
frequentieverdeling  van  een  noimale  verdeling  met  hetzelfde  gemiddelde  en  dezelfde 
standaardafwijking  is  weergegeven. 


0.50 


FORCERATIO 

Fig.  3:*  Frequendehistogiaro 

Op  basis  faiervan  kan  geconcludeerd  worden  dat  bet  twijfelacbdg  is  of  de  MOB  uit  ons  votvbeeld 
nonnaal  verdeeld  is.  Het  is  dus  nog  maar  de  vraag  of  voor  deze  dataset  een  verdelingsvrije 
methode  kan  worden  gebiuikt. 
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4.2  En^irisclie  verdelingsiunctie^ 

De  eti^iiiscbe  verdeliugsfimctie  wordt  gedefim^eTd  als 

Fq(x)  =  wr  jueer  k  van  de  in  totaal  n  waamemingen  kleiner  zijn  dan  x. 

Fg(x)  is  een  benade.ing  van  de  werkelijke  verdelingsfuncde 
F(x)  =  P(XSx)=  /'f(t)dt 

waarin  f  de  kansdichtheidsfunctie  van  X  is. 

De  volgende  figuur  geeft  de  empirische  verdelingsfunctie  van  de  MOE  uit  het  voorbeeld  in 
paragraaf  4.1,  waarbij  tevens  de  verdelingsfunctie  van  een  normale  verdeling  met  hetzelfde 
gemiddelde  en  dezelfde  variantie  is  weergegeven. 


Ook  hier  kan  geconcludeerd  worden  dat  het  twdjfelachtig  is  of  de  MOE  noraiaal  veideeld  is. 


7N0  rapport 


4.3  Statistieken 

Een  andere  manier  om  inzicht  te  krijgen  in  de  verdeling  van  een  MOE,  is  door  te  Lijken  naar  een 
aantal  statistieken.  Dit  zijn  bepaalde  gtootbeden  die  uit  de  waamemingen  worden  beiekend. 

Aantal  waamemingen 

Het  aantal  waamemingen  geeft  de  grootte  van  een  set  waamemingen  weer.  Uitgaande  van  een  set 
waamemingen  Xi..Xq,  wordt  het  aantal  waamemingen  gegeven  door  n. 

Minimum 

Het  minimum  geeft  de  laagste  waarde  uit  de  waamemingen  weer: 

Xn^n  =  (Xj) 

Maximum 

Het  maximum  geeft  de  hoogste  waarde  uit  de  waamemingen  weer; 
x^  =  mw  (Xj) 


Gemiddelde 

Het  gemiddelde  geeft  de  algemene  ligging  van  een  set  waamemingen  weer.  Veclal  wordt  onder 
het  gemiddelde  het  rekenkundig  gtmiddelde  verstaan.  De  definitie  voor  het  rekenkundig 
gemiddelde  luidt; 

x=  -  ZXj 
n  i 

Standaardafwijldng 

De  standaardafwijking  is  een  maat  voor  de  spreiding  van  de  waamemingen  rondom  het 
genuddelde.  De  definitie  voor  de  standaardafwijking  luidt: 
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Variantie 

De  variantie  is  het  kwadraat  van  de  standaardafwijking  en  wordt  gedefinieerd  door 


s2=—  I(Xi-502 

n-1  i 


De  scheefheid  is  een  maat  voor  de  symmetric  van  een  verdeling.  Een  symmetrische  verdeling 
heeft  scheefheid  0.  De  deflnitie  voor  de  schee&ieid  luidt; 


S  (Xj-^Vn 


[X  (Xj-  x)2/n]W 


De  kurtosis  gecft  de  welving  van  een  verdeling  weer.  Een  verdeling  met  een  hoge  top  en  een 
dikke  staart  heeft  een  positieve  kurtosis,  een  verdeling  met  een  lage  top  en  een  dunne  staart  heeft 
een  negatieve  kurtosis.  Een  normale  verdeling  heeft  kurtosis  0. 


kuftosis  ■  0 
kurtosis  >  0 
kurtosis  <  0 


Fig.  5:  Verdelingen  met  verschillende  kurtosis. 
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£>e  defmitie  voor  kiutosis  luidt: 


S  (Xi-x)'»/n 
i 

Y2= - 3 

[S(Xi-x)2/n]2 


Om  na  te  gaan  of  een  set  waamemingen  een  nonnale  verdeling  heeft,  is  het  zinvol  om  tc  kijken 
naar  de  scheefheid  en  de  kuitosis.  Dcze  moetcn  beiden  ongeveer  gelijk  aan  0  zijn. 


Om  te  bepalen  hoeveel  de  scheefheid  en  de  kuitosis  van  0  af  mogen  wijken,  is  1000  maal  de 
scheefheid  en  de  kuitosis  van  een  set  van  SO  waamemingen  uit  een  standaardnoimale  verdeling 
bepaald.  De  1000  waaiden  van  de  scheefheid  en  kuitosis  zijn  in  de  volgende  &equentietabel 
weergegeven: 


Frequentio  (relatief) 

Scheefheid 

Kurtosis 

-1.5  « -1.0 
-1.0  <=  -0.5 
-0.5  <=  0.0 
0.0  <=  0.5 
0.5  <=  1.0 
1.0  <=  1.5 

1.5  <=  2.0 
2.0  <=  2.5 

2.5  <.  3.0 
3.0  <=  3.5 

3.5  <=  4.0 

3  {  0.30%) 

60  (  6.00%) 
454  (  45.40%) 
427  (  42.70%) 

55  (  5.50%) 

1  (  0.10%) 

0  (  0.00%) 

0  (  0.00%) 

0  (  0.00%) 

0  (  0.00%) 

0  {  0.00%) 

14  (  1.40%) 
264  (  26.40%) 
366  (36.60%) 
220  (22.00%) 

84  (  8.40%) 

27  (  2.70%) 

16  (  1.60%) 

7  (  0.70%) 

0  (  0.00%) 

1  (  0.10%) 

1  (  0.10%) 

Dit  betekent  grofweg  dat  voor  een  set  van  SO  waamemingen  de  absolute  waaiden  van  de 
scheefheid  en  kurtosis  klelner  dan  1.0  moeten  zijn.  Is  dit  niet  het  geval,  dan  is  het  niet 
waarschijnlijk  dat  de  waamemingen  noimaal  verdeeld  zijn. 


De  waamemingen  uit  het  voorbeeld  in  paiagraaf  4.1  geven  de  volgende  statistieken: 
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FORCE_RATIO 

Aantal  waarnemingen 

50.0 

Minimum 

0.6974 

Maximum 

3.1453 

Gemiddelde 

1.4117 

Standaardafwijking 

0.5516 

Variantie 

0.3043 

Scheefheid 

1.3536 

Kurtosis 

2.0225 

De  scheefhdd  en  de  kuitosis  zijn  beiden  groter  dan  1.0,  zodat  bet  niet  waarschijnlijk  is  dat  de 
waarnemingen  uit  een  noimale  verdeling  afkomstig  zijn. 

4.4  Goodness  of  fit 

In  de  voorgaande  paiagrafen  zijn  enkele  meth^’^en  besproken  om  na  te  gaan  of  waaroemingen 
van  een  stochastische  variabele  normaal  veideeld  zijn.  Een  andere  methode  om  dit  na  te  gaan  is 
met  behulp  van  de  zogenaamde  “Goodness  of  fit"  toetsen.  Deze  toetsen  geven  weer  of  een  set 
waarnemingen  past  bij  een  bcpaalde  verdeling  (in  dit  geval  de  noimale  verdeling).  Hiervoor  kan 
een  Kolmogorov-Smimov  of  een  Chi-kwadiaat  toets  worden  gebruikt.  De  Chi-kwadraat  toets  kan 
alleen  worden  gebruikt  wanneer  de  dataset  veel  waarnemingen  bevat.  Bij  een  klein  aantal 
waarnemingen  kan  beter  de  Kolmogorov-Smimov  toets  worden  gebruikt. 

4.4.1  Kolmogorov-Smimov 

Er  zijn  twee  versies  van  de  Kolmogorov-Smimov  toets;  de  one-sample  en  de  two-sample  toets. 
De  one-sample  toets  vergelijkt  de  verdeling  van  de  waarnemingen  met  een  tbeoretiscbe  (continue) 
verdelingsfunctie.  De  two-sample  toets  vergelijkt  de  verdelingen  van  twee  sets  waarnemingen.  De 
one-sample  toets  kan  voor  goodness  of  fit  worden  gebruikt.  De  two-san^Ie  toets  komt  later  aan 
de  orde. 


In  de  one-sample  toets  wordt  de  volgende  nulbypothese  getoetst: 

Ho:F(x)  =  G(x) 
waarin: 

F(x):  de  vcrdelingsfimctie  van  de  waargenomen  variabele 
G(x):  de  verdelingsfimcde  van  de  verondersteWe  verdeling 
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£>e  toetsingsgrootheid  van  deze  toets  is 
Dnu„  =  max  ( ABS  (Fa(x)-G(x))) 

X 

waarin; 

Fn(x):  de  en^jirische  verdelingsfunctie  van  de  waargenomen  variabele. 

Deze  toetsingsgrootheid  heeft  onder  de  nulhypothese  een  bekende  verdeling.  Gegeven  een 
onbetrouwbaarheid  a  is  nu  de  krideke  waarde  R  te  bepalen  zodanig  dat 

P(D„„>R)  =  a 

Voor  n  S  35  geldt  dat 

R  >»  )  .22  /  Vn  voor  a  =  10% 

R-1.36/Vn  voor  a  =  5% 

R-»1.63/Vn  voor  a  =  1% 

Wordt  er  nu  een  toetsingsgrootheid  groter  dan  de  krideke  waarde  waargenomen,  dan  is  het  niet 
aannemelijk  (behoudens  een  onbetrouwbaarheidsmarge  a)  dat  de  toetsingsgrootheid  de 
veronderstelde  verdeling  heeft.  In  dit  geval  wordt  de  nulhypothese  dus  verworpen.  Dit  betekent 
dat  er  aangenomen  wordt  dat  de  verdelingsfuncde  van  de  waargenomen  variabele  en  de 
opgegeven  veidelingsfuncde  vcrschillen. 

Wanneer  voor  G  de  normale  verdeling  wwdt  genomen,  is  het  mogelijk  om  te  toetsen  of  de 
waamemingen  afkomstig  zijn  uit  een  normale  verdeling.  De  toets  is  echter  minder  geschikt  voor 
gevallen  waarin  de  verdeling  van  de  waamemingen  discreet  is.  Aangezicn  dit  in  het  eerder 
genoemde  voorbceld  het  geval  was,  kon  deze  toets  hier  niet  worden  gebruikt. 


TNO  rapport 


4.4.2  Chi-lcwadraat 

De  chi-kwadraat  toets  heeft  als  nulhypothese: 


Ho:F(x)  =  G(x) 


F(x);  de  verdelingsfunctie  van  de  waargenomen  variabele 
G(x):  de  verdelingsfunctie  van  de  veronderstelde  verdeling 


Bij  de  chi-kwadraat  toets  worden  de  waamenungen  opgesplitst  in  een  aantal  klassen,  aan  te 
duiden  met  n,,E2,..,Ek-  Dc  waargenomen  ftequenties  behorende  bij  deze  klassen  worden  gegeven 
door  Op  grond  van  de  in  de  nulhypothese  gegeven  verdelingsfunctie  worden  de 

firequenties  gi.gji-igtc  verwacbt.  Indien  de  nulhypothese  juist  is,  zullen  de  verschillen  tussen  de 
waargenomen  ftequenties  fj  en  de  verwachte  ftequenties  g|  klein  zijn. 

Aan  de  hand  van  deze  verschillen  kan  de  toetsingsgrootheid  worden  berekend. 


Is  de  waarde  van  deze  grootheid  klein,  dan  is  er  geen  reden  om  Hq  te  verwerpen.  Is  de  waarde  van 
groot,  dan  zal  Ho  verworpen  worden. 

Gegeven  een  onbetrouwbaaifaeid  a,  kan  de  kritieke  waarde  voor  x^  uit  de  label  van  de  chi- 
kwadraat-verdeling  worden  afgelezen.  Hierbij  moet  voor  het  aantal  vrijheidsgiaden  k-l-m  worden 
genomen,  waarin  m  bet  aantal  te  schatten  paramete's  voor  de  theoretiscbe  verdeling  is  (bij 
toctsing  tegen  een  normale  verdeling  zijn  dit  er  z).  Wordt  er  een  waarde  van  -f}  waargentmien,  die 
groter  is  dan  de  kritiei,e  waarde  uit  de  tabel,  dan  zal  de  nulhypothese  verworpen  worden. 
Voorwaarde  voor  deze  toets  is  dat  de  verwachte  ftequenties  gt  alien  groter  dan  5  zijn. 

Bij  toetsing  tegen  een  nonnale  verdeling,  zal  in  de  nulhypothese  voor  G  de  nwmale 
verdelingsfunctie  wcvden  genomen: 


v."  -  : 


TNO  rapport 


Pagina 

28 


G(x)=  J* 


1  i4((t-n)/o)2  ^ 

— i — e  dt 

<rl2n 


Het  gemiddelde  ^  en  de  standaardafwijldng  o  kimnen  uit  de  waamemingsuitlcomsten  worden 
geschat  door  het  steelqtroefgeoiiddelde  x  en  de  steekproefstandaardafwijlcmg  s.  Toetsing  van  de 
verdeling  van  X  tegen  een  normale  verdeling  met  gemiddelde  x  en  standaaidafwijking  s  komt 
overeen  met  toetsing  van  de  verdeling  van  Z  >  (X-  x)/s  tegen  een  standaard-normale  verdeling. 

Bovenstaande  procedure  werd  uitgevoerd  met  de  set  waamemingen  uit  het  voorbeeld  in  paragraaf 
4.1,  hetgeen  resulteerde  in  een  verweiping  van  de  nulhypothese.  Het  is  dus  niet  waaischijnlijk  dat 
deze  set  waamemingen  afkomstig  is  uit  een  normale  verdeling. 

4.5  Betrouwbaarbeidsintervallen 

De  ligging  van  een  reeks  waamemingen  kan  worden  gerepresenteerd  door  het 
steekproefgemiddelde  x  »  1^  I  jq.  Dit  steelqrroefgemiddelde  is  een  schatter  voor  de  verwachting 
p  van  de  verdeling  waax  de  waamemingen  uit  afkomstig  zijn.  De  bedoeling  is  nu  om  aan  de  hand 
van  de  waamemingen  X)  oen  betrouwbaaiheidsinterval  voor  p  te  constiueren.  Hierbij  ligt  het  voor 
de  hand  om  gebruik  te  maken  van  de  centrale  limietstelling: 

Gegeven  dat  de  stocbastische  variabelen  X|,..,X,  onafhankelijk  en  idendek  verdeeld  zijn  met 
verwachting  p  en  variantie  o^.  Dan  geldt  dat  wannu  n  groot  is  of  wanneer  de  Xj  noimaal 
verdeeld  zijn,  dat  x  no,  c  .aal  verdeeld  is  met  verwachting  p  en  variande 

Wanneer  er  aan  de  voorwaarden  voor  de  centrale  linuetstelling  is  voldaan,  geldt  dat 


x-p 

(Wn 


staridaard  normaal  verdeeld  is.  Hieiin  is  a  echter  nog  onbekend.  a  kan  geschat  wtnrden  door  de 
steekproefstandaardafwijldng  s.  In  dat  geval  heeft 


t  = 


x-p 

sVn 
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een  Student's  t-vetdeling  met  n-1  viijheidsgraden. 

De  lechter  kritieke  waarde  van  deze  verdeling  is  uit  de  tabel  van  de  Student's  T-verdeling  te 
halen.  Deze  wordt  aangeduid  met  to(v),  waarin  a  de  overschrijdingskans  en  v  het  aantal 
viijheidsgraden. 

Er  geldt  dus 


P(-tito(n-l)itSti4„(n-l))  =  l-a 
Dus  met  kans  1-a  geldt 


X  -  ti^(n-l)  s/Vn  n  i  X  +  tyja(n-l) 

Dit  wordt  het  1-a  betrouwbaarheidsinterval  voor  ^  genoemd. 

Bij  het  construeren  van  dit  interval  is  er  vanuit  gegaan  dat  er  aan  de  voorwaarden  van  de  centrale 
limietstelling  is  voldaan.  Dus 

nisgroot 

df  X  is  normaal  verdeeld 

Het  is  gebleken  dat  de  waaroemingen  uit  bet  voorbeeld  in  paragraaf  4.1  niet  normaal  verdeeld 
waren.  Verder  is  het  twijfelachtig  of  bet  aantal  waamemingen  van  SO  voldoende  is  om  de  centrale 
limietstelling  te  mogen  gebiuiken.  Het  is  dus  nog  maar  de  vraag  of  de  zojuist  afgeleide  foimule 
voor  een  betrouwbaarheidsinterval  voor  in  dit  geval  gebiuikt  mag  worden. 

Om  nu  toch  een  betrouwbaarheidsinterval  voor  de  verwacbting  van  deze  set  waamemingen  te 
bepalen,  wordt  gebruik  gemaakt  van  een  methode  die  geen  eisen  stelt  ten  aanzien  van  bet  aantal 
waamemingen  of  normaliteit:  de  bootstrap  methode. 


4.6 


Bootstrap  methode 


De  bootstn^  methode  biedt  de  mogelijkheid  om  voor  een  onbekende  parameter  6  met  behulp  van 
een  schatter  6  een  betrouwbaatheidsinterval  te  berekenen,  zonder  dat  daaibij  aannamen  worden 
gedaan  over  nonnaliteit  en  zonder  dat  er  eisen  worden  gesteld  aan  bet  aantal  waamemingen. 
Hierbij  woidt  gebruik  gemaakt  van  Mon*  *  Carlo  simulatie. 

Uitgaande  van  een  steekproef  met  n  waamoningen,  worden  er  uit  deze  steekproef  n  trekldngen 
met  teruglegging  gedaan.  De  aldus  verkregen  set  van  waamemingen  wordt  een  bootstn^ 
steekproef  genoemd.  Deze  procedure  wordt  1000  maal  herhaald.  Vervolgens  wordt  voor  elke 
bootstrap  steelqrroef  de  bootstn^  schatter  8*  bepaald  (het  sterreQ'e  geeft  aan  dat  de  schatter 
betrekldng  heeft  op  een  bootstrap  steekproef). 

De  verkregen  bootstrap  schattingen  kunnen  in  een  frequendehistogram  worden  uitgezet.  Dit  geeft 
een  beeld  van  de  verdeling  van  S.  Om  nu  een  1  -  a  betrouwbaaifaeidsinterval  voor  6  te  berekenen, 
wordt  gebmik  gemaakt  van  de  percentile  methode  (zie  [Efron  2]). 

Zij  G(s)  de  verdelingsfimcde  van  , 

G(s)  =  P(0*^s) 

Dan  is  [G'>('/6a},  G'>(l'Via)]  een  benadering  voor  bet  1  •  a  betrouwbaaifaeidsinterval  voor  8.  Om 
de  percentielen  G'KVia)  en  G'>  (1-Via)  te  bepalen,  wordt  G  benadeid  door  de  empirische 
verdelingsfuncde  Giooq.  Om  het  Via  percendel  van  Gjooq  te  bepalen,  worden  de  waargenomen 
bootstrap  schattingen  geordend.  De  geordende  bootstrap  schattingen  worden  ^2*>  ••  >  ^looo* 
genoemd.  Wanneer  voor  de  k*  bootstrap  schatdng  8^*  geldt  dat 

Giooo(  ®k*)  ^  <  Giooo(  ®k+l*) 

A 

dan  is  8^*  bet  Via  percendel  van  Gjgoo.  Op  dezelfde  manier  kan  bet  1  •  Via  percendel  van  Giggo 
worden  bepaald.  Het  aldus  gevonden  interval  geeft  een  benadering  van  het  1  -  a 
betrouwbaadieidsinterval  voor  8. 

Een  betrouwbaaifaeidsinterval  voor  de  verwacbdng  |i  lf«n  nu  verkregen  worden  door  bet 
bovenstaande  toe  te  passen  met  als  schatter  bet  gemiddelde. 
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4.7  Datasets  vagelijken 

Er  zijn  verschillende  statisdsche  methoden  om  sets  waamemingen  met  elkaar  te  vergelijken.  Deze 
methoden  woiden  onderscheiden  in  verdelingsafliankelijke  en  verdelingsvrije  methoden. 

4.7.1  Verdelingsafhankelijke  methoden 

Enkele  verdelingsafhankelijke  methoden  zijn  de  Student's  t-toets,  de  toets  van  Satterthwaite  en  de 
F'toets  voor  vaiianties.  Bij  deze  methoden  wordt  er  aangenomen  dat  de  waamemingen  normaal 
verdeeld  zijn.  De  Student's  t-toets  en  de  toets  van  Satterthwaite  toetsen  of  twee  sets 
waamemingen  dezelfde  verwachting  hebben,  de  F-toets  voor  varianties  toetst  of  twee  sets 
waamemingen  dezelfde  verwachting  hebben. 

4.7. 1 . 1  Student's  t-toets 

De  Student's  t-toets  lean  wotden  gebruikt  om  te  toetsen  of  er  versclul  in  de  verwachtingen  van 
twee  sets  waamemingen  is.  De  toets  heeft  als  voorwaarde  dat  de  waamemingen  uit  beide  sets 
normaal  verdeeld  zijn,  en  dat  de  varianties  van  beide  verdelingen  geUjk  zijn. 

Gegeven  zijn  2  sets  waamemingen: 

X„X2...,X„ 

yi.y2.”.ya 


Wanneer  Xi..Xm  afkomstig  zijn  uit  een  verdeling  met  verwachting  en  variantie  en  yi..yQ  uit 
een  verdeling  met  verwachting  en  variantie  dan  is  de  nulhypothese  voor  de  t-toets: 


De  toeisingsgrootbeid  is 


t* 


*-y 


waarin  s  dc  woctcl  is  van  het  gewogen  gcmiddcldc  van  dc  sicekprocfvariantics  cn  : 
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s 


(m-1)  S)2  +  (n-1) 
m  +  n-2 


De  toetsingsgrootheid  heeft  onder  de  nulhypothese  een  Student's  t-verdeling  met  m+n-2 
vrijheidsgraden.  Uit  de  tabel  van  de  Student's  t-verdeling  zijn  nu  de  kiitieke  waaiden  en  R2  te 
halen  zodanig  dat 


P(t<R,)*=V4a 
P(t>R2)  =  Wa 

Wordt  er  een  toetsingsgrootheid  kleiner  dan  R|  of  groter  dan  Rj  waargenomen,  dan  wordt  de 
nulhypothese  vcrworpen. 

Ook  is  het  mogelijk  om  een  1-a  betrouwbaarheidsinterval  voor  |i)-^2  te  constnieren ; 


x-y  -  t«a(m+n-2)  s  A/  -  +  -  <  Hj-Hj  <  x-  y  +  ti5a(m+n-2)  s  \J  -  +  - 
V  m  n  V  m  n 

Ligt  de  waarde  0  binnen  dit  interval,  dan  is  er  geen  reden  om  de  nulhypothese  te  verwerpen. 

4.7. 1 .2  Toets  van  Satterthwaite 

De  t-toet$  mag  alleen  toegepast  worden,  wanneer  er  aan  de  voorwaarde  van  gelijke  varianties  is 
voldaan.  Is  dit  niet  bet  geval,  dan  kan  er  een  aangepaste  versie  van  de  t-toets  worden  gebiuikt  die 
deze  voorwaarde  niet  heeft.  Deze  toets  wordt  de  toets  van  Satterthwaite  genoemd. 
Toetsingsgrootheid  voor  deze  toets  is 


**y 

.  So 


waarin 
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Onder  de  nulhypothese  benadcrt  de  verdeling  van  dcze  toetsingsgrootheid  de  t-verdcling  met  f 
vrijheidsgradcn,  waarbij  f  wordt  gcgeven  door 


m-1  n-1 


Op  dezelfde  manier  als  bij  de  t-toets  kunnen  er  nu  kritieke  waaiden  worden  bepaald  en 
betrouwbaarheidsintervallen  worden  berekend. 

4.7. 1 .3  F-toets  voor  varianties 

De  F-toets  voor  varianties  kan  gebruikt  worden  om  te  toetsen  of  twee  sets  waamemingen  dezelfde 
variantie  hebben.  Wanneer  de  eerste  set  bestaat  uit  ni  waamemingen  uit  een  noimale  verdeling 
met  verwachting  |Xi  en  variantie  Oj^,  en  de  tweede  set  uit  n  waamemingen  uit  een  noimale 
verdeling  met  verwachting  Hj  en  variantie  dan  is  de  nulhypothese  voor  dezc  toets: 


Ho :  (ii^  X  O2* 


De  toetsingsgrootheid  voor  deze  toets  is 


waaiin  Sj^  Tt  $2^  (is  dit  nict  het  geval,  verander  dan  de  nummeiing  van  beide  sets).  Onder  de 
nulhypothese  heeft  F  een  F-verdeling  met  m-1  en  n-1  vrijheidsgraden.  De  kritieke  waarde  voor  F 
is  nu  uit  de  tabel  van  de  F-verdeling  af  te  lezen,  en  bij  overschrijding  van  deze  waarde  wmrdt  de 
nulhypothese  verworpen. 

4.7.2  Verdelingsvrije  methoden 

Enkele  verdelingsvrije  methoden  zijn  de  Kolmogorov-Smimov  two-sample  toets,  de  toets  van 
Wilcoxon  en  een  variant  van  de  bootstrap  methode.  Deze  methoden  doen  geen  aaimamen  over  de 
verdeling  van  de  waamemingen.  De  Kolmogorov-Snumov  two-sample  toets  gaat  rut  of  twee  sets 
waamemingen  uit  dezelfde  verdeling  afkomstig  zijn.  De  toets  van  Wilcoxon  en  de  variant  van  de 
bootstrap  methode  toetsen  of  twee  sets  waamemingen  Jezelfde  verwachting  hebbot. 


4.7.2.1  Kolmogorov-Smimov  (two-san^le) 

De  Kolmogorov-Smimov  two-sample  toets  kan  gebruilct  woiden  om  te  toetsen  of  het 
waarscfaijnlijk  is  dat  twee  sets  waamemingen  uit  dezelfde  verdeling  afkomstig  zijn.  Hieibij  is  het 
niet  noodzakelijk  dat  beide  sets  evenveel  waamemingen  hebben.  De  two-sample  toets  heeft  als 
nulhypothese 


Ho;F(x)  =  G(x) 


woarin 

F(x):  de  verdelingsfhnctie  van  de  eerste  set  waamemingen 
G(x);  de  verdelingsfunctie  van  de  tweede  set  waamemingen 

De  nulhypothese  is  dus  dat  beide  datasets  dezelfde  verdelingsfunctie  hebben. 


De  toetsingsgrootheid  van  deze  toets  is 


Da„«  max  (abs  (Fn(x)  -  Go,(x))) 
% 


waarin  Fg(x)  de  empirische  verdelingsfunctie  van  de  eerste  datase..  en  G^Cx)  de  empirische 
verdelingsfunctie  van  de  tweede  dataset. 

Deze  toetsingsgrootheid  heeft  onder  de  nulhypothese  een  bekende  verdeling.  Gegeven  een 
onbetrouwbaarheid  a  is  nu  de  kritieke  waarde  R  te  bepalen  waarvoor  geld* .  U 

P(D^>R)»o 

Wanneer  er  een  toetsingsgrootheid  groter  dan  de  kritieke  waarde  wordt  waargenomen,  wordt  de 
nulhypothese  verwotpen.  Het  is  dan  niet  aannemelijk  (behoudens  een  cmbetrouwbaatheidsmarge 
a)  dat  beide  sets  waamemingen  uit  dezelfde  verdeling  afkomstig  zijn. 

In  plaats  van  te  toetsen  of  de  verdelingifuncties  gelijk  zijn,  is  het  ook  nsogelijk  om  te  toetsen  of  de 
verdelingsfuncties  niet  meer  dan  een  bepaalde  tolerantie  6  van  elkaar  veischillen.  Deze  toets  heeft 
als  nulhypothese 


Ho:F(x-5)5G(x)SF(x+5) 


De  naar  rechts  verschoven  fimctie  F(x-S)  wordt  vericregen  door  bij  de  oorsprookelijke 
waamemingen  S  op  te  tellen,  de  naar  links  verschoven  funcde  F(x+5)  wordt  vedcregen  door  van 
de  oorspronkelijke  waamemingen  5  af  te  trekken.  De  volgende  figuur  geeft  een  voorbeeld  van 
twee  verdelingsiimcties  die  niet  voldoen  aan  de  bypothese  F(x)  =  G(x),  maar  voor  een  bepaalde  S 
wel  voldoen  aan  Hg: 


X 

Fig.  6:  Voorbeeld 

Eenvoudig  is  aan  te  tonen  dat  de  nulhypothese  Hg  overeenkomt  met 


Hg:G(x-5)SF(x)SG(x+6) 

zodat  het  niet  uitmaakt  welke  verdelingsfunctie  F(x)  wordt  genoemd,  en  welke  verdelingsfimctie 
G(x). 

Is  de  tolerande  6  gelijk  aan  0,  dan  geeft  dit  de  al  ecrder  beschrcven  two-san^le  toets. 


De  toetsingsgrootheid  voor  de  toets  met  tolerande  is 
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Doa*  =  niax(D+,D-) 

waarin;  D+=  max  (Fn(x-5)  -  Gn,(x)) 

X 

D-=  max(G„(x)-F„(x+S)) 


VooT  deze  toetsingsgrootheid  lean  weer  de  krideke  waaide  worden  bepaald.  Wordt  er  een 
toetsingsgrootheid  groter  dan  de  kritidee  waarde  waargenomen,  dan  wordt  de  nulhypothese 
verwoipen.  Dit  houdt  in  dat  het  niet  aannemelijk  is  dat  de  functie  G(x)  tussen  F(x>$)  en  F(x-f5)  in 
Ugt. 

Opmerldng:  voor  6  ongelijk  aan  0  is  niet  bewezen  dat  dezelfde  verdeling  heeft  als  voor  S 
gelijk  aan  0.  Dit  is  echter  wel  aannemelijk. 

AH  22  Toets  van  Wilcoxon 

De  toets  van  Wilcoxon  kan  gebruikt  worden  om  te  toetsen  of  er  verschil  is  tussen  de  ligging  van 
twee  sets  waamemingen.  Stel  de  eerste  set  waamemingen  is  afkorosdg  uit  een  verdeling  met 
verwachdng  jij,  en  de  tweede  set  is  afkomstig  uit  een  verdeling  met  verwachdng  Ma-  ^^n  is  de 
nulhypothese  voor  deze  toets 


Ho:^l  =  H2 

Voor  het  bepalen  van  de  toetsingsgrootheid  van  deze  toets  worden  de  waamemingen  van  beide 
sets  in  een  geordende  tij  gezet  (beide  sets  door  elkaar).  Vervolgens  worden  aan  de  elenaenten  van 
deze  rij  rangnummers  Qi.CJz-iQomi  toegekend.  De  toetsingsgrootheid  is  nu  de  som  van  de 
rangnummers  van  de  waamemingen  uit  de  eerste  dataset: 


W=  2  Q, 

id 

Is  de  waarde  van  deze  toetsingsgrootheid  klein,  dan  zijn  de  waamemingen  uit  de  eerste  dataset 
kleiner  dan  de  waamemingen  uit  de  tweede  dataset.  Is  de  waarde  van  deze  toetsingsgrootheid 
groot,  dan  zijn  de  waamemingen  uit  de  eerste  set  groter  dan  die  uit  de  tweede.  In  beide  gevallen 
zal  de  nuliiypothese  verwotpen  worden.  Er  zal  dus  tweezijdig  getoetst  moeten  worden. 
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De  toetsingsgrootheid  heeft  onder  de  nulhypothese  een  bekende  verdeling.  Gegeven  een 
onbetrouwbaarheid  a  kunnen  nu  de  kritieke  waaiden  Ri  en  R2  worden  bepaald  zodanig  dat 

P(W<R,)  =  Via 
P(W>R2)  =  V4a 

Zijn  m  en  n  groot  genoeg,  dan  is  W  nonnaal  verdeeld  met  verwachting 

1 

Hw  ®  ~  (m+n)  ♦  (m+n+1) 

4 

en  variantie 


SB  —  m  *  n  *  (m+n+1) 

12 

De  kritieke  waarden  Ri  en  R2  worden  dan  gegeven  door 

Ri  *  ^  *  Ow 

R2  =  Hw  +  Zyta  *  Ow 

Voor  kleine  m  en  n  moeten  deze  waarden  in  een  tabel  worden  opgezocht.  Wordt  er  een 
toetsingsgrootheid  kleiner  dan  R|  of  groter  dan  R2  waargenomen,  dan  wordt  de  nulhypothese 
verworpen. 

4.7.2.3  Bootstrap  variant 

De  bootstrap  methode  (zie  4.6)  kan  wrsden  gebiuikt  om  een  betrouwbaarbeidsinterval  voor  het 
verschil  in  de  vtawachtingen  van  twee  verdelingen  te  berekeneu.  Hieibij  wordt  als  schatter  het 
verschil  tussen  de  steekproefgemiddelden  genomen. 
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Gegeven  zijn  2  sets  waamemingen 


Ui,U2,..,Un 

Vl.V2...,Vn 

Wanneer  Uj..Ug,  afkomstig  zijn  uit  een  verdeling  met  verwachting  en  Vi..Vq  uit  een  verdeling 
met  verwachting  |X2,  dan  is 


-  2uj  -  -  Svj 

m  t  n  J 

een  zuivere  schatter  voor  Nu  kan  er  met  behulp  van  de  bootstn^  methode  een  benadering 
voor  bet  1-a  betrouwbaarheidsinterval  voor  jli-jl2  worden  beiekend.  De  bootstrap  steekproeven 
die  hicrbij  gebruikt  worden,  bevatten  m  trekkingen  met  teruglegging  uit  U|,..,Ubi  en  n  trekldngen 
met  teruglegging  uit  Vi,..,Vn. 

Wanneer  0  binnen  het  betrouwbaarheidsinterval  ligt,  is  er  geen  reden  om  aan  te  nemen  dat  pj  en 
P2  verschillen.  Ligt  0  niet  binnen  dit  interval,  dan  zal  de  hypothese  dat  en  P2  zijn  moeten 
worden  verwotpen.  In  dit  geval  zijn  er  twee  mogelijkheden; 

1 .  Het  interval  ligt  geheel  boven  0  hetgeen  implicecrt  dat  p j  groter  is  dan  H2 

2.  Het  interval  ligt  geheel  beneden  0  hetgeen  impliceert  dat  Pi  kleiner  is  dan  P2 
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5  DEO^rIWIKKELDEPROGRAMMATUIJR 

In  het  van  het  berekenen  van  de  MOE's  en  het  uitvoeren  van  stadstische  toetsen  zdjn  drie 
conqjuterprogranuna's  gescbreven.  Deze  zijn  geprogrammeerd  in  Ada  en  diaaien  onder  VAX- 
VMS.  Deze  drie  prognumna's  zijn; 

-  WRTTE.UNITS 

Dit  programma  haalt  de  informade  uit  de  uitvoer  van  het  Force  Stiuchire  Model  (FSM)  die 
nodig  is  om  de  MOE's  te  berekenen.  De  uitvoer  van  dit  progranuna  is  een  tekstfile  met 
gegevens  over  eenheden. 

-  CALCULA're_MOE 

Dit  programma  berekent  waarden  voot  een  aantal  MOE's  aan  de  hand  van  gegevens  over 
eenheden.  De  uitvoer  is  een  matrix  van  getallen  waarin  per  simuladerun  de  «aarden  voor 
een  aantal  MOE's  staan. 

-  DEC_MOE_VIS 

Dit  is  een  interacdef  programma  dat,  aan  de  hand  van  (6n  of  meerdere  matrices  met 
gegevens,  stadsdsche  analyses  uitvoert. 

Later  in  dit  hoofdstuk  zuUen  deze  drie  programma's  uitgebreid  worden  beschreven.  Daarbij 
worden  voorbeelden  van  uitvoer  en  invoer  van  de  programma's  weergegeven  door 
Vette  letters  geven  een  keuzemogelijkheid  weer  van  het  interaedeve  programma 
DEC_MOE_VIS.  Vanwege  afspraken  die  gemaakt  zijn  over  de  uitvoer  van  programmatuur  die  in 
het  leader  van  FSM  worden  geschreven  is  deze  uitvoer  en  de  interactive  interface  in  het  Engels. 

De  uitvoer  van  het  ene  programma  dient  als  invoer  voor  een  ander  progranuna.  Dit  is  gedaan  om 
de  aihankelijkheid  van  het  FSM  model  kwijt  te  raken.  Ook  op  andere  manieren  kan  invoer 
gegcnereeid  worden  voor  de  programma's  CALCULATE_MOE  en  DEC_MOE_VIS.  Figuur  7 
geeft  schemadsch  aan  hoe  de  invoer  van  de  verschillende  programma's  verkregen  kan  worden; 
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Fig.  7:  Scbematisch  overzicht  invoer  verschillende  progranuna's 


Het  programma  WRrnB_UNrrS  verkrijgt  .»ijn  invoer  uit  bet  gevecbtssimulatieinodel  FSM.  Uit  dc 
uitvoer  van  WRITE_UNrrS  worden  met  het  programma  CALCULATE_MOE  verschUlende 
MOE's  berekend.  Dit  programma  kan  ook  MOE's  voor  uitkomsten  van  andere 
gevechtssimulatiemodellen  berekenen,  mits  deze  uitkpmsten,  eventueel  met  behulp  van  een  filter, 
in  het  juiste  formaat  staan.  Met  behulp  van  het  programma  DEC_MOE_VIS  kunnen  de  berekende 
MOE's  vervolgens  worden  geanalyseerd.  Dit  programma  kan  ook  gebruikt  worden  om  de 
uitkomsten  van  een  willekeurig  stochastisch  simulatiemodel  te  analyseren.  Hiertoe  moet  aan  de 
uitkomsten  van  het  simulatiemodel  een  waardering  worden  toegekend. 
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5.1  Het  computcrprogramma  WRITE_UNITS 

Het  programma  WRITE.UNITS  haalt  uit  de  uitvoer  van  het  FSM  model  de  gegevens  waarmee 
f  verschillende  MOE's  berekend  kunnen  worden.  Voordat  het  programma  gedraaid  wordt,  moct  er 

in  een  aparte  file  staan  op  welk  moment  van  de  simulatie  de  gegevens  weggeschreven  moeten 

I 

worden.  Hiervoor  zijn  verschillende  mogelijkheden; 

Alij  de  toestand  op  een  bepaalde  tijdsdp  in  de  simulatie  van  belang  is  (bijvoorheeld  de 
toestand  op  10: 15:00  uur)  dan  komt  in  de  file  te  staan: 

TIME*:  10:15:00 

Als  de  toestand  van  belang  is  bij  het  bereiken  van  een  bepaalde  'marker',  dan  komt  in  de  file 
te  staan: 

MARKER- 111  DOOR 

Deze  markers  worden  door  het  FSM  model  in  de  uitvoerfile  geschrcven.  Tijdens  een 
simulatie  kan  door  middel  van  condities  in  het  scenario  gezorgd  worden  voor  het  leggen  van 
markers.  Zo  kan  er  bijvoorbeeld  een  marker  gelegd  worden  op  het  moment  dat  de  eenheid 
1 1 1  door  een  mijnenveld  is  gebroken. 

Als  de  toestand  van  belang  is  waarop  de  sttrkte  van  een  bepaalde  eenheid  of  een  combinatie 
van  eenheden  beneden  een  op  te  geven  drempel  komt.  Bijvoorbeeld  de  regel 
BLUE  30% :  4NT  +  dllM  +  4N  AT  <  4 

geeft  aan  dat  de  toestand  wordt  weggeschreven  op  het  moment  dat  het  totaal  aantal 
voertuigen  van  4NT,  4  NM  en  4NAT  kleiner  dan  of  gelijk  is  aan  4.  BLUE  30%  geeft  een 
naain  aan  zodat  later  kan  worden  teruggevonden  dat  de  toe.stand  om  deze  reden  is 
weggeschreven. 

In  de  file  kunnen  meerdere  regels  staan  waarop  de  toestand  moet  worden  weggeschreven.  De 
toestand  wordt  echter  maar  d<n  keer  weggeschreven  en  wel  zo  vrocg  mogelijk.  Als  bijvoorbeeld 
in  de  file  staat: 

RED40%:11<16 

BLUE  30% :  4NT  +  4NM  +  4N  AT  <  4 

TIME  -  10:20:00 


Dan  wordt  de  toestand  weggeschi  }ven 
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als  het  aantal  voertuigen  van  de  eenheid  1 1  kleiner  dan  of  gelijk  is  aan  16,  of 

als  het  gezamenlijk  aantal  voertuigen  van  4NT,  4NM  en  4NAT  kleiner  dan  of  gelijk  is  aan  4, 

of 

als  het  tijdstip  10:20:00  is  bereikt,  of 

als  de  siniulatie  gestopt  is  voordat  €6a  van  de  eerder  genoemde  voorwaarden  bereikt  is. 


Van  elke  simulatienin  tvordt  bepaald  op  welk  tijdstip  en  om  welke  reden  de  toestand  wordt 
weggeschreven.  Met  het  progranuna  DEC_MOE_VIS  kan  een  ovendcht  gegeven  worden  hoe 
vaak  de  toestand  om  een  bepaalde  reden  is  weggeschreven  (de  dumpreden). 

De  set  van  dumpredenen  staat  in  een  gewone  tekstfile.  Er  mogen  ook  commentaarregels  en  lege 
regels  voorkomen  in  deze  file.  In  de  file  RSM_USER^FILES  moet  de  logical  DUMP_SET 
aangeven  in  welke  file  de  set  van  dumpredenen  staat.  De  parameters  FIRSTJIUN.ID  en 
LASTJllJN_ID  uit  de  file  FSM_USER_PARAMETERS  geven  aan  van  welke  runs  de  toestand 
moet  worden  weggeschreven.  Meer  informatie  over  de  files  FSM_USER_FILES  en 
FSM_USER_PARAMETERS  staat  in  bijlage  B. 

De  uitvoer  van  het  progranuna  WRITE.IJNITS  bestaat  uit  een  aantal  tekstfiles.  Voor  de 
beginsituatie  en  voor  elke  run  wordt  een  aparte  uitvoerfile  geaetierd.  De  naam  van  de  files  wordt 
gemaakt  aan  de  hand  van  de  logical  UNTT.DATA,  het  versienummer  van  FSM  en  het 
runiK'TJner.  In  de  eerste  regel  staat  informatie  over  de  gehele  run  en  in  de  andere  regels  staat 
informatie  over  de  eenheden.  De  regels  met  aan  het  begin  van  de  regel  zijn  commentaarregels 
en  zijn  geen  informatie.  Ze  zijn  bedoeld  om  de  uitvoer  leesbaarder  te  maken. 

De  eerste  regels  uitvoer  van  dit  programma  zouden  bijvoorbeeld  kunnen  zijn; 


—  I^UN^IP 

DONP  MUSOM 

OOKP.TXMX 

101 

RED  401 

IOjTtiOO 

— OHXTJHMa 

COtOR 

wirjgxft 

wi_or_vtH 

B08XT20M 

4KAn 

SLOS 

TOWl” 

o" 

0  0 

4KATI 

TONS 

0 

0  0 

BLOK 

TOW2 

0 

0  0 

BLOC 

TOW! 

0 

0  0 

4 

BLOK 

XKBC 

0 

0  0 

BLOK 

wa 

1 

SMIlO  SBB7640 

iMC 

BLOC 

WQ 

1 

5»iB70  Bt$7410 

Uit  de  eerste  regel  blijkt  dat  dit  de  uitvoer  is  van  run  101,  dat  de  toestand  weg  is  geschreven 
omdat  aan  de  voorwaarde  RED  40%  werd  voldaan.  Dit  gebeurde  om  10:17:00  uur. 
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5.2  Het  programnu  CALCULATE_MOE 

Het  programma  CALCULATE_MOE  boekent  aan  de  hand  van  de  gegevens  van  de  eenbeden  de 
waarden  van  een  aantal  MOE's.  De  invoer  is  ecn  tekstfile.  De  uitvoer  van  WRITE_UNITS  kan  als 
invoer  gebruikt  worden,  maar  eventueel  kan  ook  op  andeie  wijze  een  dergelijke  tekstfile 
aangemaakt  worden.  De  gegevens  van  de  eenheden  moeten  dan  echter  wel  op  dezelfde  manier  in 
een  file  staan  als  de  gegevens  in  de  uitvoeifile  van  WRrrE_lINlTS. 

De  eisen  die  het  programma  CALCULATE_MOE  oplegt  aan  de  invoer  zijn: 

0  Op  de  eerste  regel  staan  de  RUN_E>,  de  naam  van  de  DUMP _REASON  en  het  tijdstip 
waarop  de  toestand  gedumpt  is. 

De  RUN_ID  is  een  identificatie  van  de  run  waaniit  deze  gegevens  komen.  Dit  is  een  geheel 
getal. 

De  naam  van  de  DUMP_REASON  is  een  naam  die  aangeeft  om  welke  reden  deze  toestand 
gedumpt  is.  Deze  naam  is  10  karakters  lang. 

Het  tijdstip  van  dumpen  wordt  aangegeven  in  het  formaat  UREN:M1NUTEN:SECX>NDEN. 

0  Op  elke  andere  regel  staat  informatie  over  6€a  eenheid.  Deze  informatie  bestaat  uit: 

De  naam  van  de  eenheid  (maximaal  10  karakters,  geen  spades) 

De  kleur  van  de  eenheid  (BLUE  of  RED) 

Het  type  van  de  eenheid  (maximaal  10  karakters,  geen  spades) 

•  Het  aantal  voertuigen  van  deze  eenheid  op  het  djdsdp  waarop  t  e  toestand  is  gedumpt 
(geheel  getal) 

De  positie  van  de  eenheid.  Deze  wwdt  aangegeven  door  een  x>cootdinaat  cn  een  y- 
coordinaat.  Beiden  zijn  gehele  getallen  en  zijn  gegeven  in  meters.  Bijvootbeeld  UTM- 
codrdinateu  kunnen  worden  gebruikt.  De  poside  (0, 0)  wordt  gegeven  bij  eenheden  die 
geen  voertuigen  mecr  hebben. 

Lege  ibgels  cn  commentaar  (aangegeven  door  --)  worden  oveigeslagen  bij  het  inlezen. 

In  de  file  FSM_USER_FILES  moet  aangegeven  worden  welke  dies  gebruikt  moeten  worden  als 
invoer  en  als  uitvoer.  De  logical  UNIT_DATA  geeft  de  luuun  van  de  invoodles  aan.  Deze  naam 
wordt  automadsch  aangevuld  met  een  vetsienuramer  en  een  runjd.  De  logical  MOE_DATA 
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geeft  aan  hoe  de  uitvoerfile  heet.  In  de  file  FSM_USER_PARAMETERS  geven  de  parameters 
FIRST_RUN_ID  en  LAST_RUN_ID  aan  welke  invoerfiles  gcbruikt  worden. 

De  MOE's  die  momenteel  berekend  worden  zijn  (zie  hoofdstuk  3  voor  een  be^.chrijving  van  de 
MOE’s): 


-  FORCE.BLUE 

-  LOSS.RED 

-  FORCEJIATIO 

-  LOSS.RATIO 

-  Wm.FACTOR 

-  obst_defenc:e 

De  tijdsduur  tussen  aanvang  van  het  scenario  en  moment  van  dumpen.  Op  rich  is  dit  niet 
echt  een  MOE  (een  tijdsduur  van  20  is  niet  altijd  beter  of  slechter  dan  een  tijdsduur  van  IS), 
maar  in  combinatie  met  de  dumpreden  kan  dit  wel  een  MOE  zijn. 

De  overige  MOE's  worden  niet  berekend  omdat  dit  niet  gemakkelijk  en  eenduidig  kan.  Voor 
specifieke  gevallen  ka.n  het  programma  CALCULATE_MOE  aangepast  worden  zodanig  dat  deze 
MOE's  wel  uitgerekend  kunnen  worden. 


uitvoer  van  CALCULATE_MOE  wordt  geschreven  in  een  tekstfile.  De  file  krijgt  de  naam 
zoals  die  in  de  FSM_USER_FILES  door  de  logical  MOE_DATA  wordt  gedefinieerd.  De  eerste 
en  laatste  regels  zien  er  bijvoorbeeld  als  voIgt  uit; 


0  —  rirvt 
7  —  tast  »o«_ind«x 

lOl  —  rirat  Tun_ld 
150  —  laat  XMn^id 


DUMP  MASON 

3.00000C>00 

Z.OOOOOX^OO 

1. 00000X^00 

2.00000K«00 

2»OOOOOK«00 

2.00000X400 


roxcx^itox 

5.(2SOOK-Ol 

6,175001-01 

3.7S000X-01 

6.35000X-01 

5.62500K-01 

7.50000K-01 


tOSS^XXD 

S.I0744X-01 

5.II744X-01 

5.9I463K-01 

5.S9744K-01 

5.t9744X-01 

5.i9744X-01 


rOXCKMtlO 

1.347t»K«00 

1.  ••710X400 

•.61539K-01 

1.57265X400 

1.547tfX400 

2.35097X400 


LOfSJATXO 

1.37109X400 

1.67570X400 

•.12500X-01 

1.52344X400 

1.37109X400 

1.02013X400 


WIW^rxCTOR 

1.35954X400 

1.70140X400 

0.37019X-01 

1.54004X400 

1.35954X400 

2.09355X400 


OMTJOXrXM 

7.71406X-01 

0.136G7X-01 

3.69921X-01 

7.50532X-01 

6.696IOK-01 

0.92701X-01 


DOTATION 

1.70000X401 

1.65000X401 

2.30000X401 

7.50000X400 

1.65000X401 

7.50000X400 


2.OC0OOX4OO  6.97900X-01  5.09744S-01  1.00710X400  1.67570X400  1.70140X400  0.96036X-01  1.60000X401 
2.00000X400  5.00000C-01  5.09744S-01  1.17949X400  1.31075X400  1.19912X400  6.03202X-01  1.70000X401 
XLOX  301  CKtlf»«r«MQO  1 
RXO  404  —  Dua^saaaoa  2 


De  getallen  op  de  eerste  vier  regels  geven  bet  fonnaat  aan  van  de  nutria  met  MOE*waarden.  In 
elke  rij  staan  de  resultaten  van  simulatierun.  EUce  kolom  van  de  matiix  geeft  de  waarden  van 
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66n  MOE.  De  kolommen  hebben  dtels.  Dit  zijn  de  namen  van  de  betreffende  MOE's.  In  de  eeiste 
kolom  van  de  matrix  staan  getallen  die  aangeven  om  welke  leden  de  toestand  is  weggeschreven. 
Onderaan  staan  de  dtels  voor  de  dumpredenen. 

5.3  Het  programma  DEC_MOE_VIS 

Het  programma  DEC_MOE_VIS  is  een  interacdef  programma  waatmee  een  aantal  stadsdsche 
grootheden  uitgerekend  kan  worden  en/of  gradsche  afbeeldingen  gemaakt  kunnen  worden  van  de 
waarden  van  de  MOE's.  Eerst  zal  beschreven  worden  hoe  met  het  programma  gewerkt  moet 
worden.  Dit  zal  gedaan  worden  aan  de  hand  van  een  aantal  voorbeelden.  Daama  zal  een 
beschrijving  van  de  benodigde  invoer  volgen.  Het  programma  maakt  gebniik  van  de  IMSL 
bibliotheek  voor  een  aantal  stadsdsche  methoden  en  maakt  gebniik  van  MONGO  voor  het 
tekenen  van  de  dguren. 

5.3.1  Handleiding  programma  DEC_MOE_VIS 

Na  het  opstarten  verschijnen  er  menu's  met  een  aantal  keuzemogelijkheden.  Door  middel  van  de 
pijltjestoetsen  omhoog  en  omlaag  wordt  de  gewenste  mogelijkheid  gekozen.  Het  pijltje  naar  links 
betekent  terug  naar  het  vorige  menu  en  met  return  wordt  een  bcpaalde  keuze  uitgevoerd. 

Na  het  opstarten  verschijnt  het  volgende  menu  in  beeld; 

Force  structure  Model  V12 


ilUi  -  Vit.  I  l«-fk  »  111..-  J  /-AUi.- 1  '  ' )  “ 


Fig.  8:  Hoofdroenu 

AUereerst  zal  voor  keuzemogelijkheid  7  (Setup)  gekozen  rooeten  worden.  Hietmee  wmdt 
ingesteld  waar  de  tii^^oer  naar  toegaat.  De  ui'.voer  is  te  verdelen  in  twee  klassem 
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De  grafische  uit''oer.  Hierbij  kan  gdcozen  worden  uit  een  aantal  mogelijke  typen  appaiaten. 
Door  middel  van  het  menu  kan  het  gewenste  type  uitvoeiapparaat  gekozen  worden.  Daama 
zal  het  device  aangegeven  moeten  worden.  Dit  kan  een  terminal  zijn  of  een  filenaam 
(afhankelijk  van  het  gekozen  type).  Kiest  men  voor  NONE  dan  wordt  er  geen  grafische 
uitvoer  gegenerecrd. 

De  tekst  uitvoer.  Kiest  men  in  het  setup  menu  voor  te;it  dan  moet  men  daama  de  filenaam 
invoeren  van  de  tekstfile  waar  de  uitvoer  naar  toe  geschreven  moet  worden.  Bij  het  vragen 
naar  de  filenaam  verschijnt  als  standaard  antwoord  "NONE".  Laat  men  dit  staan  als 
filenaam,  dan  zal  er  geen  tekst  uitvoer  gegenereerd  worden. 

Na  het  instellen  van  de  uitvoer  parameters  kunnen  de  stadstische  opties  gekozen  worden. 

5.3.2  Voorbeelden  van  het  werken  met  DEC_MOE_VIS 

In  deze  paragraaf  wordt  met  behulp  van  een  aantal  voorbeelden  een  oveizicht  gegeven  van  de 
mogeUjkheden  van  DECJvIOE_VlS.  Als  iets  vet  gediukt  is  dan  geeft  dit  een  keuzemogelijkheid 
weer  van  het  programma  DEC 3I0E_V1S  en  als  iets  ki»in  gedrukt  is  dan  is  dit  een  voorbeeld  van 
de  uitvoer  zoals  die  in  de  tekstfile  komt  te  staan  die  met  Setup  en  Text  is  aangegeven.  De 
voorbeelden  die  in  deze  paragtaaf  worden  beschreven  zijn: 

Voorbeeld  1 :  Frequentiehistogram  met  gegevens  uit  66n  dataset. 

Voorbeeld  2;  Frequentiehistogram  met  twee  datasets. 

Voorbeeld  3:  Empirische  verdelingsiunctie  van  twee  datasets 
Voorbeeld  4:  Statistieken  van  66n  dataset 
Voorbeeld  5:  Vergelijken  dataset  met  uoimaie  verdeling 
Voorbeeld  6:  Betrouwbaarheidsinterval  voordeverwacbting 
Voorbeeld  7:  Het  vergelijken  van  twee  normaal  verdeelde  datasets. 

Voorbeeld  8:  Het  vergelijken  van  twee  datasets  met  behulp  van  de  Kolmogorov>Smimov  toets 
Voorbeeld  9:  Her  vergelijken  van  twee  datasets  met  behulp  van  de  toets  van  Wilcoxon 
Voorbeeld  10:  Het  vergelijken  van  twee  datasets  met  behulp  van  de  Bootstrap  roetbode 
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Voorbeeld  1:  Frequendehistogram  met  gegevens  uit  ^  dataset. 

De  gegevens  zijn  gesplitst  naar  hun  dunqtreden.  Op  deze  manier  kan  zichtbaar  gemaakt  worden  of 
er  veel  verschil  is  tussen  de  verdelingen  van  de  uitkomsten  bij  de  diverse  dumpredenen. 

Frequencytable 

De  gebruika  wil  een  frequentiehistogram. 

1  set  with  dunip_reasons 

De  gebmiker  wil  een  firequentiehistogram  waarbij  de  gegevens  uit  66n  set  komen.  De 
gegevens  zijn  opgesplitst  naar  een  dumpreden. 
dataset  1 :  $dirl;moe_data 

De  gegevens  van  de  betieffende  MOE  komen  uit  deze  file, 
dataset  with  dump_reasons :  $dirl:moe_data 

De  gegevens  betreffende  de  indeUng  van  de  waamemingen  naar  dumpreden  staan  in  deze 
die  (dit  hoeft  niet  noodzakelijkerwijs  dezelfde  file  te  zijn  als  de  file  met  de  MOE's). 

Moe  choice :  forceratio 

De  gebmiker  wil  een  frequentieverdeling  van  de  MOE  forceratio. 


De  figuur  bestaat  uit  twee  delen: 

1  In  het  onderste  gedeelte  staan  de  uitk(»nsten  per  run.  Op  de  X-as  staat  de  uitkomst  voor  de 
gekozen  MOE  en  op  de  Y-as  een  index  voor  het  nummer  van  de  run.  Een  steireQ'e  geeft  aan 
dat  de  uitkomst  uit  een  set  waamemingen  komt  met  alien  de  eerste  dunqneden  en  een 
driehoekje  duidt  op  een  uitkomst  uit  een  set  met  de  tweede  dunqireden.  Bijvoorfoeeld  de 
uitkomst  van  runnummer  101  is  1.33  en  heeft  als  dumpreden  RED  40%. 

2  De  bovenste  helft  geeft  een  reladeve  frequentiehistogram.  Deze  geeft  bijvoorbeeld  aan  dat 
18%  van  de  waamemingen  uit  de  set  met  de  eerste  dumpreden  komt  en  tussen  0.S  en  1.0  ligt 
en  dat  4%  van  de  waamemingen  uit  de  set  met  de  tweede  dumpreden  komt  en  tussen  0.5  en 
1.0  ligt. 

De  bijbehorende  uitvoer  in  de  tekstfile  ziet  er  als  volgt  uit: 

$01RllH0E  DATA 


rOKCERXTXO 

Fraquaney 

(zalativa) 

BLUE  20% 

BIO  40% 

0,5  <-  1.0  : 

9 

(  IS.OOt) 

2  (  4.00%) 

1,0  <-  1.5  t 

0 

(  0.00%) 

23  (  46.00%) 

1.5  <-  2.0  t 

0 

(  0.00%) 

11  <  22.00%) 

2.0  <•  2.5  i 

0 

(  0.00%) 

3  (  6.00%) 

2,5  <«  5.0  t 

0 

(  0.00%) 

0  (  0.00%) 

3,0  <•  3.5  X 

0 

(  V.00%) 

2  (  4.00%) 

Total 

9 

(  15.00%) 

41  (  52.00%) 

Uit  de  tabel  valt  af  te  lezen  dat  er  3  waamemingen  zijn  in  de  set  met  de  tweede  dumpreden  en  met 
een  waarde  tussen  2.0  en  2.S,  ofwel  6%  van  alle  waamemingen. 

Het  aantal  klassen  en  de  klassegrenzen  worden  door  het  programma  bepaald  aan  de  hand  van  de 
minimum-  en  maximumwaarde  van  de  waamemingen. 

Voorbceld  2:  Frequentiehistogram  met  twee  datasets. 

Frequency  table 

De  gebruiker  wil  een  fc.quentiehistogram. 

2  sets 

De  gebniiker  wil  de  frequendes  van  twee  datasets  met  elkaar  vergelijken. 

Dataset  1:  $dirl:moe_data 

Dit  is  de  file  met  gegevens  van  de  eerste  data  set 
Dataset  2:  $dii2:moe_dau 

Dit  is  de  file  met  gegevens  van  de  tweede  dau  set 
Moe  choice :  forceiado 

De  gebruiker  wil  een  frequentieverdeling  van  de  MOE  forcerado. 


0.50 

0.^5 

0.^0 

0.35 
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De  figuur  bestaat  uit  twee  delen: 

1  In  het  onderste  gedeelte  staan  de  uitkomsten  per  run.  Op  de  X-as  staat  de  uitkomst  voor  de 
gekozen  MOE  en  op  de  Y-as  tec  xodex  voor  het  nummer  van  de  run.  Ecn  steire^e  geefl  aan 
dat  de  uitkomst  uit  de  eerste  data&le  komt  en  een  driehoelge  duidt  op  ecn  uitkomst  uit  de 
tweede  datafile.  Bijvoorbeeld  de  uitkomst  van  runnummer  101  uit  de  eerste  dataset  is  1.3S 
en  uit  de  tweede  dataset  is  2.35. 

2  De  bovenste  helft  geeft  een  relatieve  frequentiehistogram.  Deze  geefl  bijvoorbeeld  aan  dat 
46%  van  de  waamemingen  uit  de  eerste  dataset  tussen  1.0  en  1.5  ligt  en  dat  33%  van  de 
waamemingen  uit  de  tweede  dataset  tussen  1.0  en  1.5  ligt. 

De  bijbehorende  uitvoer  in  de  tekstfile  ziet  er  als  volgt  uit: 

rORC£RATXO  rr«qu«ncy  {r«l*tlv*) 


S«t  1 

s.t 

2 

0.5  <- 

1.0  : 

11 

<  22.004) 

11  ( 

11.004) 

1.0  <- 

1.5  j 

23 

(  46.001) 

33  ( 

33.004) 

1.5  <• 

2.0  : 

11 

<  22.004) 

41  ( 

41.004) 

2.0  <- 

2.5  1 

3 

(  6.004} 

9  ( 

9.004) 

2.5  <- 

3.0  : 

0 

<  0.004) 

0  ( 

0.004) 

3.0  <- 

3.5  : 

2 

<  4.004) 

4  ( 

4.004) 

3.5  <- 

4.0  : 

0 

(  0.004) 

0  ( 

0.004) 

4.0  <- 

4.5  t 

0 

<  0.004) 

0  ( 

0.004) 

4.5  <- 

5.0  i 

0 

(  0.004) 

2  ( 

2.004) 

1  •  $DXR1:H0£_DATA 
S«t  2  •  $0XR2iH0S*DATA 

Uit  bovenstaande  label  blijkt  dat  23  waamemingen  uit  de  eerste  set  (46%  van  de  waamemingen 
uit  de  eerste  set)  een  waarde  hebben  tussen  1.0  en  1.5. 

Voorbeeld  3:  Empirische  .i.verdelingsfunctie;  van  twee  datasets 

Cumulative  distribution 
2  sets 

De  gebruiker  wil  de  verdelingsfuncties  van  twee  datasets  met  elkaar  vergelijken. 

Dataset  1;  $diil;moe_data 

Dit  is  de  file  met  gegevens  van  de  eerste  data  set. 

Dataset  2:  $dir2:moe_data 

Dit  is  de  file  met  gegevens  van  de  tweede  data  set. 

Moe  choice ;  forceratio 

De  gebruiker  wil  de  verdelingsfuncties  van  de  MOE  forceratio. 


0.9 

0.8 

0.7 

0.6 
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De  flguur  bestaat  uit  twee  delen: 

1  In  het  onderste  gedeelte  staan  de  uitkomsten  per  nm.  Op  de  X-as  staat  de  uitkomst  voor  de 
gekozen  MOE  en  op  de  Y-as  een  index  voor  het  nrnnmer  van  de  run. 

2  De  bovenste  helft  geeft  de  verdelingsfuncties.  Deze  geeft  bijvoorbeeld  aan  dat  10%  van  de 
waamemingen  uit  de  eerste  dataset  kleiner  of  gelijk  is  aan  0.7  en  dat  10%  van  de 
waameming.n  uit  de  tweede  dataset  kleiner  of  gelijk  is  aan  0.8. 

Voorbeeld  4:  Statistieken  van  6dn  dataset 


Statistics 

De  gebruiker  wil  statistieken  van  een  dataset. 

Dataset  1:  $dirl:moe_data 

Dit  is  de  naani  van  de  ^e  met  de  gegevens 
Dataset  with  dump_reasons :  $dirl:inoe_data 

De  gegevens  betreffende  de  indeling  van  de  waamemingen  naor  dumpreden  staan  in  deze 
file  (dit  hoeft  niet  noodzakelijkerwijs  dezelfde  file  te  zijn  als  de  file  met  de  MOE's). 

Moe  choice:  forceratio 

De  gebruiker  wil  de  sutistieken  van  de  verdeling  van  de  MOE  forceratio  berekenen. 


De  uitvoer  in  de  tekstfile  ziet  er  als  volgt  uit: 


roRccm^Tio 

Ml  oht«rvtttion» 

M«An 

1.4117 

Varl4nc«  : 

0.3043 

dt4nd4id  d«vi4tion  : 

0.5314 

8k«vn«s»  i 

1.353< 

Rurtosl*  1 

2.0225 

Klniaua  : 

0.6974 

KaxIbub 

3.1453 

Nx.  ob*«rv4tion* 

50.0 

BLUE  301 

N4&n  i 

0,7977 

Varitnc*  : 

0.0047 

Standard  d«vi«tion  i 

0.06S4 

Sk«wn«s« 

0.1025 

Kuttoala  x 

-1.1967 

KlnlBua  X 

0.6974 

HaxlmuB  i 

0.9026 

Nr.  obsarv'Ationa  x 

9.0 
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Mean 

1.S464 

VariancA 

0.2684 

Standard  daviation 

o.siei 

Skawnasa 

1.S893 

Kurtosis 

2.3636 

HiniiBun 

0.9436 

Maximum 

3. 1453 

Mr.  obsarvatlona 

41.0 

De  eerste  tabel  geeft  de  statistieken  van  alle  waamemingen  en  de  volgende  twee  tabellen  geven 
de  statistieken  van  de  waarnemingen  gesplitst  naar  de  dumpreden.  Skewness  is  de  Engelse 
benaming  voor  scheefheid. 


Voorbeeld  5:  Vergelijken  dataset  met  nonnale  verdeling 
Goodness  of  fit 

De  gebruiker  wil  verdeling  vergelijken  met  normale  verdeling. 

Chi'Square 

De  gebruiker  wil  een  goodness  of  fit  test  doen  met  behulp  van  de  chi-kwadraat  toets. 

Dataset  1:  $dirl:moe_data 

Dit  is  de  naam  van  de  file  met  de  gegevens 
Level  of  confidence:  95.0 

Het  betrouwbaarheidsniveau  1  -  a  is  95% 

Moe  choice:  forceratio 

De  gebruiker  wil  de  verdeling  van  de  MOE  forceratio  vergelijken  met  een  normale 
verdeling. 


De  uitvoer  in  de  tekstfile  ziet  er  als  volgt  uit: 


Chi'squar*  . . . . 

0At«««t  1 

$DIIU  }HOC_DATA 
Mo«  ChoiCA 


rORCERATlO 

t«v«l  of  confidAnc*  ;  95.0000 
Chl-aquACAd  statistic  :  17.1996 
P'valuA  ;  0.0006 
or  :  3 , 0000 
Th#  hypothssis  will  ba  rAjsetad 


j 


\ 

I 

I 


«  - 
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Dc  P-value  geeft  aan  hoe  groot  de  leans  is  dat,  als  de  nulhypothese  waar  is,  cr  cen 
toetsingsgrootheid  gevonden  wordt  die  groter  of  gelijk  is  aan  de  uitgerekende  chi-kwadraat 
toetsingsgrootheid.  Als  deze  P-value  kleiner  is  dan  a,  dan  wordt  de  nulhypothese  verworpen.  DF 
is  het  aantal  viijheidsgraden. 

Voorbeeld  6:  Betrc’iwbaaiheidsinterval  voor  de  verwachting 


Confidence  intervals 

De  gebruiker  wil  een  betrouwbaarheidsinterval  voor  de  verwachting  berekenen. 
Dataset  1;  $dirl;tnoe_data 

Dit  is  de  naam  van  de  gewenste  dataset 
Level  of  confidence;  95.0 

Het  betrouwbaarheidsniveau  1  -  a  is  95% 

Moe  choice;  forceratio 

De  gebruiker  wil  het  betrouwbaarheidsinterval  van  dc  MOE  forceratio  berekenen. 


De  uitvoer  in  de  tekstfile  ziet  er  als  volgt  uit: 


Con£i<S«nc«  lnt«rv«lf 
1 

^DlfU  :H0E_0ATA 
Ho«  Choic« 
rCRCSRATlO 


Ccnfl4«nc«  int«rv4la  for  tb«  population  iima  yaluo 


Lovol  of  confidonco  : 
Nofm!  nathod  : 
Bootstrap  standard  asthod  t 
Bootstrap  paresntila  asthod  ; 
Bias'Ccrzsetsd  psrcantils  aathod: 


95.0000 
1.2549  * 
1.2449  * 
1,2727  - 
1.2915  - 


1.5494 

1.5595 

1.5433 

1.5701 


Wanneer  het  betrouwbaarheidsinterval  volgens  de  normale  methode  wordt  beiekend,  krijgt  nten 
als  antwoord  dat  de  verwachting  met  een  waarschijnlijkbeid  van  95%  ligt  tussen  1.2549  en 
1.5684.  Dit  is  echter  met  de  aanname  dat  de  MOE  normaal  verdeeld  is.  De  andere 
betrouwbaarheidsintervallen  zijn  met  behulp  van  de  bootstrap  methode  verkregen. 


V 

A 

t 

I 

t 
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Voorbeeld  7:  Het  vergelijken  van  twee  nonnaal  verdeelde  datasets. 


Compare  datasets 

De  gebniiker  wil  datasets  vergelijken. 

Student's  t 

De  gebruiker  wil  de  Student's  t*toets  of  daaimee  samenhangende  toetsen  gebruiken. 
Dataset  1:  $dirl;moe_data 

Dit  is  de  naam  van  de  eerste  dataset 
Dataset  2;  $dir2;moe_data 

Dit  is  de  naam  van  de  tweede  dataset 
Level  of  confidence;  95.0 

Het  betrouwbaarhcidsniveau  1  -  a  is  95% 

Moe  choice;  forceratio 

De  gebruiker  wil  de  verwachtingen  van  de  MOE  forceratio  vergelijken. 


De  uitvoer  in  de  tekstiile  ziet  er  als  volgt  uit; 


Cosptit*  4«t4««ta 

8tud«nt'«  t  . . 

1 

$D1RI  :H0E_DATA 
2 

$01R2  tH0E_CA7A 
Mo«  Choice 
rORCCFATIO 

of  confid«nc«  t  95.0000 

Confidanc*  interval  for  tha  dlffaranca  in  saana  (first  alnua  aacond)  ( 
ASSUM2NO  EQUAL  VARIANCES  : 

-0.4552  -  -0,0154 

Tha  hypothaala  (aaana  ara  aqual)  vlXl  ha  rajactad  I 

KOT  ASSUKXNG  EQUAL  VARIANCES  ; 

-0.4406  -  -0.0300 

Tha  hypothaala  (maana  ara  aqvAl)  vlll  ba  rajactad  ! 

Kypothaaa:  varlancaa  ara  aqual 
Laval  of  confldanca  :  95.0000 

r-valua  {  1.5325 

Da9raaa  of  fraadoa  t  49«  99 

Prob.  of  a  lar^ar  P  :  0.0992 

Tha  hypothaala  will  ba  rajactad  1 


Deze  procedure  tnag  allecn  gebruikt  worden  als  aangenoroen  mag  worden  dat  beide  sets 
waaraemingen  uit  een  normale  verdeling  komen.  Uit  dit  voorbeeld  blijkt  dat  met  een  kans  van 
95%  het  verschil  tussen  de  twee  verwachtingen  ligt  tusscn  -0.4552  en  -0.0154  als  ook 
aangenomen  mochl  worden  dat  de  varianties  aan  elkaar  gelijk  zijn.  Mag  de  aanname  van  geUjke 
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vaiianties  niet  gedaan  worden  dan  ir  bet  betrouwbaadteidsinterval  in  dit  voorbeeld  -0.44Q6  tot  -  t 

0.0300  (berekend  met  de  Satterthwaite  procedure). 

Ook  wordt  een  toetsingsgrootheid  (F-value)  gegeven,  waarmee  de  hypothese  getoetst  kan  worden 
dat  de  twee  varianties  aan  elkaar  gelijk  zijn  (F-toets).  De  kans  op  een  grotere  F  moet  vergeleken 
worden  met  de  onbetrouwbaarfaeid  a.  Is  de  kans  kleina  dan  a,  dan  dient  de  nulhypothese 
verworpen  te  worden.  Het  is  af  te  raden  om  deze  toets  te  gebruiken  om  te  bepalen  welk  van  de 
twee  gegeven  betiouwbaarheidsintervallen  gebtuikt  wordt.  Bij  twijfel  over  gelijkheid  van  de 

! 

variandes  moet  Satterthwaite's  procedure  worden  ^rbniikt. 

Voorbeeld  8:  Het  vergelijken  van  twee  datasets  met  behulp  van  de  Kolmogorov-Smimov  toets 

Compare  datasets 

De  gebruiker  wil  datasets  vergelijken. 

Kolmogorov-Smirnov 

De  gebruiker  wil  de  Kolmogorov-Smimov  toets  gebruiken. 

Dataset  1:  $dirl;tnoe_data 

Dit  is  de  naam  van  de  eerste  dataset 
Dataset  2;  $dir2;moe_data 

Dit  is  de  naam  van  de  tweede  dataset 
Level  of  confidence;  95.0 

Het  betrouwbaarheidsniveau  1  •  a  is  95% 

Tolerated  difference:  0.0 
De  tolerande  is  0.0 
Moe  choice:  forccrado 

De  gebruiker  wil  de  verdelingen  van  de  MOE  forcerado  vergelijken. 

De  uitvoer  in  de  tekstfile  ziet  er  als  volgt  uit: 

KoIa.09orov*SalrnoY . 

1 

fDXIUJNOC^DATA 
D«t4S«t  2 
$OXlC:HOe_DATA 
M04  Cholca 
rORCXAATZO 

of  ccnfldaaco  :  9S.OOOO  ' 

diiXtssonc*  i  v<vvvv 

0+  *  0.2400 

t>-  *  0.0000  *  j 

P-V4l«*  I  0.0299 

Tb4  bypoth«4i»  will  b«  f  • 
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D(-  P-value  geeft  aan  hoe  groot  de  leans  is  dat,  als  de  nulhypothese  waar  is,  er  een 
toetiingsgrootheid  gevonden  wordt  die  groter  of  gelijk  is  aan  de  gevonden  toetsingsgrootheid 
MAX  (D+,  D-).  Als  deze  P-value  kleiner  is  dan  a,  dan  wordt  de  nulhypothese  verwoipen. 

In  dit  voorbeeld  wordt  de  nulhypothese  dat  de  verdelingen  van  de  MOE’s  gelijk  zijn  verwotpen. 

Voorbeeld  9:  Het  vergelijken  van  twee  datasets  met  behulp  van  de  toets  van  Wilcoxon 


Compare  datasets 

De  gebruiker  wil  datasets  vergelijken. 

Wilcoxon 

De  gebruiker  wil  de  toets  van  Wilcoxon  gebruiken. 

Dataset  1;  $dirl;moc_data 

Dit  is  de  naam  van  de  eerste  dataset 
Dataset  2;  $du2:moe_data 

Dit  is  de  naam  van  de  tweede  dataset 
Level  of  confidence:  95.0 

Het  betrouwbaarheidsniveau  1  -  a  is  95% 

Moe  choice:  Lreeratio 

De  gebruiker  wil  de  verwachtingen  van  de  MOE  forceratio  vergelijken. 


De  uitvoer  in  de  tekstfile  ziet  er  als  voIgt  uit: 


Coaf>«r« 

Wilcoxon . . 

1 

9DX1U  :NOK_OATA 
0«t«x«t  2 

|02A2iM0K_DA7A 
Moo  Cholco 
rOACTMTlO 

Lovol  of  confldoACo  t  95.0000 
Wilcoxon  tost  •totlstle  t  3134.5000 


Tbo  hypotboslo  will  bo  rojoctod  t 


De  P-value  geeft  aan  hoe  groot  de  kans  is  dat,  wanneer  de  nulhypothese  waar  is,  er  een 
toetsingsgrootheid  gevonden  wordt  die  groter  of  gelijk  is  aan  de  uitgerekende  Wilcoxon 
toetsingsgrootheid.  Als  deze  P-value  kleiner  is  dan  a,  dan  wordt  de  nulhypothese  verwotpen. 
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Voorbeeld  10:  Het  vergelijken  van  twee  datasets  met  behulp  van  de  Bootstrap  methode 

Compare  datasets 

De  gebruiker  wil  datasets  vergelijken. 

Bootstrap 

De  gebruiker  wil  de  bootstrap  methode  gebiuiken. 

Dataset  1:  $dirl;moe_data 

Dit  is  de  naam  van  de  eerste  dataset 
Dataset  2;  $dir2;moe_data 

Dit  is  de  naam  van  de  tweede  dataset 
Level  of  confidence;  95.0 

Het  betrouwbaarheidsniveau  1  -  a  is  95% 

Moe  choice:  forceratio 

De  gebruiker  wil  de  verwachtingen  van  de  MOE  forceratio  vergeUjken. 

De  uitvoer  in  de  tekstfile  ziet  er  als  volgt  uit: 


BootiCrap . . . 

Oataaat  X 

90XR1  iHO£_DA7A 
Datatat  2 

9DXR2  >H0£_DA7A 
Moa  Cholca 
rORCCBATZO 

Laval  of  confidanca  s  9&.0000 

Confldanca  Intarval  foe  tha  diffaranca  In  saana  (flrat  minua  aacond)  : 
-0.4313  -  -0.0344 

Tha  hypothaala  (aaana  ara  aquaU  vill  ba  rajactad  I 


Uit  dit  voorbeeld  blijkt  dat  met  een  kans  van  95%  het  verschil  tussen  de  twee  verwachtingen  ligt 
tussen  -0.4313  en  -0.0364.  Omdat  0.0  niet  in  dit  interval  ligt  is  wordt  de  nulhypothese  (beide 
verwachtingen  zijn  aan  elkaar  gelijk)  verworpen. 
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5.3.3  De  benodigde  invoer 

De  invoer  van  het  programma  DEC_M0E_V1S  bestaat  uit  een  matrix  met  daarin  de  waarden  voor 
de  MOE's.  Op  eUce  regel  (rij)  staan  de  uitkomsten  van  66n  run.  In  elke  kolom  staan  de  uitkomsten 
van  66n  MOE.  In  de  file  met  deze  matrix  worden  eerst  enkele  getaUen  verwacht  die  de  grootte 
van  de  matrix  weergeven.  Verder  moeten  de  namen  van  de  MOE's  in  de  IBle  staan  voor  de  matrix. 
Een  vootbeeld  van  de  eerste  en  laatste  regels  van  een  invoerfile  ziet  er  als  volgt  uit: 


0  "*•  rir«t  no«^lnd«x 
7  — ‘  Last  noa^indax 

101  —  rirat  run^ld 

150  Last  run  id 


DUKPRCXSON 

2.00000E400 

2.00000C400 

1.00000C400 

2.00000E400 

2.00000E400 

2.OC000E4O0 


rORCE_BLOE 

S.C2S00E*01 

6.87S00E-01 

3.75000E>01 

e.2S000E-01 

S.<2500E«01 

7.S0000E-01 


LOSS^REO 

S.$9744£-01 

5.89744E-01 

S.38462E-01 

S.I9744E-01 

5.89744E*01 

5.8d?44E<01 


rORCERRTlO 

1.34799E400 

1.68718E400 

8.61S39E-01 

1.57285E400 

2.347ft9C400 

2.3S8m400 


LOS8_RMZO 

1.37To9E400 

1.67S78S400 

8.12500E-01 

1.52344E400 

1.37109E400 

1.82813S400 


HlH^rACTOR 

1.35954S400 

1.78148E400 

8.37019E-01 

1.548041400 

1.359S4E400 

2.09355E400 


OEST^DEriN 

7.7140IE-01 

8.136071-01 

3.635211-01 

7.58532E-01 

6.69680C-01 

8.92701E-01 


DURATION 

1.70000E401 

1.65000E401 

2.30000E401 

7.500001400 

1.6S000E401 

7.50000B400 


2.00000E400  6.87S00E-01  S.89744E-01 
2.00000E400  S.OOOOOE-01  5.897441-01 
BLUE  301  —  Dump-raaaon  1 
RED  40^  -•  Oufflp-raaacn  2 


1.88718E400  1.67578E400  1.781481400  8.968361-01 
1.179491400  1.218751400  1.199121400  6.832821*01 


1.600001401 

1.700001401 


De  eerste  vier  regels  geven  de  grootte  van  deze  matrix  wecr.  Ei  zijn  8  kolommen  in  de  matrix  (0 
tot  en  met  7)  en  SO  rijen  (101  tot  en  met  ISO).  In  kolom  0  staat  een  getal  dat  aangeeft  om  welke 
dumpreden  de  toestand  is  weggeschteven.  Als  er  verschillende  dumpredenen  gebruikt  zijn  moet 
dit  altijd  ^la  kolom  nul  aangegeven  worden.  De  nummcrs  van  de  rijen  (101  tot  en  met  150)  horen 
bij  de  i  un-id's  van  de  verschillende  runs. 

De  vijfde  regel  geeft  de  namen  van  de  kolommen,  dus  van  de  MOE's.  Een  naam  bestaat  uit 
maximaal  10  karakters  en  mag  geen  spades  bevatten. 

Na  de  matrix  met  reele  getallen  kutmen  enkele  regels  volgen  met  op  elke  regel  een  naam  van  een 
dumpreden. 

De  uitvoer  van  het  programma  CALCULATE_MOE  kan  als  invoer  dienen  voor  dit  programma. 

In  FSM_USER_FILES  moet  aangegeven  worden  welke  files  als  invoer  gebruikt  moeten  worden. 
Dit  wordt  gedaan  door  de  logicals  MOE_DATA_l,  MOE_DATA_2  enz.  In  bijlage  B  staat  meer 
over  bet  gebruik  van  FSM_USER_FILES. 
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6  MOGELDKHEDENVOORVERDERONDERZOEKENAANBEVELINGEN 

Nu  er  de  mogelijkheden  zijn  om  op  eenvoudige  wijze  maten  van  effecdviteit  (MOE's)  uit  te 
rekenen  en  te  analyseren,  zdjii  er  ook  uitgebreidere  statistische  ondeizoeken  mogelijk. 

De  mogelijkheden  tot  ve.'der  ondeizoek  die  in  dit  hoofdstuk  behandeld  worden  zdjn; 
Gevoeligheidsanalyse 
Metamodellen 


Gevoeligheidsanalyse 

Gevoeligheidsanalyses  worden  uitgevoerd  om  te  bepalen  of  het  model  gevoelig  is  voor  kleine 
veranderingen  in  de  invoer.  Stel  voor  een  invociparameter  een  interval  vast  waartussen  deze 
parameter  kan  schommelen.  Neem  nu  de  twee  uitersten  en  ga  daaimee  scenario's  doorrekenen.  De 
andere  invoerparameters  worden  constant  gehouden.  De  gevoeligheid  kan  worden  bepaald  door 
de  uitkomsten  met  elkaar  te  vergelijken.  Zie  paragraaf  4.7  (Datasets  vergelijken). 

Het  kan  nodig  zijn  om  enkele  sets  van  invoerparameters  door  te  rekenen,  omdat  in  66a  set  de 
uitkomst  van  het  model  ongevoelig  kan  zijn  voor  deze  invociparameter  terwijl  in  een  andere  set 
wel  gevoeligheid  optreedt.  Als  bij  66n  van  de  sets  de  uitkomsten  gevoelig  zijn  voor  deze 
invoeiparameter  dan  is  het  model  gevoelig  voor  deze  parameter. 

Metamodellen 

Bij  metamodelleiing  worden  relades  (F)  geschat  tussen  een  set  invoerparameters  (xj, ..  en  de 
verwachtingswaarde  van  d<n  of  meerdere  maten  van  effectiviteit  (MOE). 


MOE  =  F(x„..,ig 


Zo  kan  bijvoorbeeld  een  relatie  geschat  worden  tussen  de  initiele  gevechtskrachtverhouding  en 
verwachtingswaarde  van  de  terreinwinst  van  de  vijand.  Deze  lelaties  kunnen  globale  inzichten 
verschaffen  in  de  gevoeligheid  voor  veranderingen  in  de  invoerparameters: 


MOEj  =  F  (x,. .. ,  Xi  +  5, ...  Xo) 
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Ook  kunnen  de  relades  gebniikt  worden  om  globale  schatdngen  te  maken  van 
verwachtingswaarden  van  MOE's  met  waaiden  van  invoeiparameters  waaivoor  geen  simulaties 
zijn  gedraaid: 


MOE'  =  F(x,’,..,x„') 

Om  de  schatdng  van  de  relatie  F  te  kunnen  maken  worden  de  waarden  van  de  invoerparameters 
gevarieerd.  Uit  theoretisch  oogpunt  wordt  verondersteld  dat  de  simulatie-uitkomst  a£bankelijk  is 
van  deze  invoeiparameters. 

Voor  elke  set  invoerwaarden  worden  een  aantal  simulaties  gedraaid  en  aan  de  hand  van  de 
simulatie-uitkomsten  kan  een  schatting  gemaakt  worden  van  de  verwachting  van  de  MOE. 
Doordat  er  meerdere  sets  zijn  met  waarden  van  de  invoerparameters  en  bijbeborende  schatdngen 
van  de  verwachting  kan  de  relatie  geschat  worden  tussen  de  invoeiparameters  en  de  verwachting. 

Om  dit  sooit  analyses  te  kunnen  doen  is  het  nodig  dat  er  goede  inzichten  zijn  in  de  theoretische 
samenhang  tussen  invoeiparameters  en  de  gekozen  MOE. 

Verder  is  een  opzet  nodig  voor  het  uitgebreid  beschrijven  van  de  invoer.  Als  er  een  goede  opzct  is 
en  bij  elk  scenario  worden  tevens  gegevens  genoteerd  over  de  gevonden  MOECs)  dan  kan  na 
verloop  van  tijd  onderzoek  gestart  worden  naar  verbanden  tussen  waarden  van  invoerparameters 
en  de  uitkomst  van  de  MOE. 
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HET  GEBRUIV  VAN  FSM_USER_FILES  EN  FSM_USER_PARAMETERS 

In  deze  uijlage  wordt  beschreven  hoe  de  files  FSWLUSERJE’ILES  en 
FSM_USER_PARAMETERS  emit  zien.  Deze  files  zijn  om  parameters  en  filenamen  in  tc  stellen 
die  afhankelijk  zijn,  of  kunnen  zijn  van  de  omgeving.  Bijvoorbeeld  de  namen  van  invoer-  en 
uitvoeifiles  of  de  indices  van  de  runs  die  gedraaid  worden.  Het  gebniik  van  deze  files  impliceert 
dat  ook  een  aantal  andere  files  aanwezig  moeten  zijn.  Hieronder  staan  alle  benodigde  files 
weergegeven  en  is  aangegeven  wat  er  tninimaal  in  moet  staan. 


DefileFSM  SYSTEM  FILES. 

In  deze  file  staan  verwijzingen  naar  de  andere  files.  De  file  kan  er  als  volgt  uitzien: 


FSM  SYSTEM  PARS  FILE 

FSM  USER^PARS  FILE 

ERROR^FILE  FILE 

EML.FILE  FILE 

LOG.TERM  TERMINAL 


QFSM.SYSTEM  PARAMETERS 

QFSM.USER.PARAMETERS 

DMOE.ERRORS 

QERRSLOG 

nl: 


EML_FILE  geeft  de  logical  voor  de  logfile  en  LOG_TERM  voor  de  log_tcnninal.  In  deze  file  en 
op  deze  tenninal  kunnen  door  de  programma's  meldingen  en  waarscbuwingen  worden 
weggeschreven  die  het  mogelijk  maken  in  de  gaten  te  houden  wat  het  programma  doet.  De 
programma's  die  in  dit  rapport  beschreven  worden  maken  nauwelijks  gebruik  van  deze 
mogelijkheid.  Het  is  dan  ook  aan  te  raden  om  de  logical  LOG.TERM  te  definieren  als  nl:,  dit 
betekent  dat  er  geen  log_terminal  wordt  geopend. 


DefileFSM  SYSTEM  PARAMETERS. 

In  deze  files  staan  parameters  die  in  alle  progranuna's  noodzakelijk  zijn.  De  file  kan  er  als  volgt 
uit  zien: 


RUN_NAME  STRINGJTYPE  MOE 

START_MODE  START.MODES  SCRATCH 


FSP_LOG_LOGFILE  BOOLEAN_TYPE  FALSE 

FSP_LCX}_TERMINAL  BOOLEAN  TYPE  FALSE 

FSP_AC_LOG_LOGFILE  BOOLEAnItYPE  FALSE 

FSP_AC_LOG_TERMINAL  BOOLEAN  TYPE  FALSE 
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De  laatste  vier  panuneter^  geven  aan  of  er  meldingen  over  het  gebruik  van  de  parameters 
weggeschreven  worden  in  de  logfile  of  op  de  lo£_terminal. 


Defde  PARAMETER  DATABASE. 

In  deze  file  staat  een  ovendcht  van  de  aangeroepen  parameters  en  files  door  de  verschillende 
programma's. 


DefileMOE  ERRORS.DIRECT 

In  deze  file  staan  de  foutmeldingen  die  in  het  programma  DEC_MOE_VIS  gegenereerd  Inmnen 
worden. 


PefileFSM  USER  FE.ES. 

Deze  files  geven  aan  welke  files  als  invoer  en  uitvoer  voor  de  verschillende  programma's  worden 
gebruikt.  Onderstaand  voorbeeld  geeft  aan  hoe  deze  file  emit  zou  kunnen  zien.  Achter  elke 
logical  staat  aangegeven  welk  programma  deze  logical  nodig  heeft.  Als  een  ander  programma 
gedraaid  wordt  kan  deze  logical  weggelaten  worden. 


DBWI  FILE 

FILE 

$DATDIR.XXXX_ 

-WRITE  UNITS 

UDT_FILE 

FILE 

DUSER_DEFINED_TYPES 

-WRTTE.UNITS 

PDT.FDLE 

FILE 

SDATDDLPRE  DEFINED  TYPES 

-  WRITE.UNITS 

ANALYSIS_FILE  FILE 

□ANALYSIS 

-WRITE  UNITS 

MUF_FILE 

FILE 

□UNITS 

-WRITE.UNrrS 

RIF_FILE 

FILE 

□RIF 

-WRITE  UNITS 

DUMP_CRrrERIAFILE 

□DUMP.TXT 

-WRITE_UNITS 

UNTT.DATA 

FILE 

□UNIT.DATA 

-  WRITE.UNrrS 
&  CALCULATE_MOE 

GRID  TERRAIN 

FILE 

$DATDIR:GRID  TERRAIN 

-  CALCULATE  MOE 

MOE.DATA 

FILE 

□MOE.DATA 

-  CALCULATE_MOE 

MOE  DATA  1 

FILE 

$DIR1:M0E_DATA 

-DECJdOE  VIS 

MOE_DATA^2 

FILE 

$DIR2:MOE_DATA 

-  DECJklOE.VIS 

$DATDIR,  SDIRl 

en  $DIR2  zijn  in  dit  voorbeeld  logicals  die 

in  de  VMS  omgeving  zijn 

gedtfiniecrd.  Het  grid-tenein  wordt  allecn  gebmikt  om  de  MOE  OBST_DEFEN  (hindemis 


verdedigingskracht)  uit  te  rekenen. 
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DefileFSM  USER  PARAMETERS 

In  deze  file  staan  enkele  parameters  die  in  eokele  van  de  programma's  gebruikt  worden.  In 
onderstaand  voorbeeld  wordt  aangegeven  in  welk  programma  de  parameters  gebruikt  worden. 


RUN.ID  INTEGER.TYPE  100 

FIRST.ID  INTEGER_TYPE  101 


LAST_ID 


INTEGER_TYPE  150 


_  ALLE 

~WRrrE_UNITS 
&  CALCULATE.MOE 
-  WRITE_UNITS 
&  CALCULATE  MOE 


START_TIME 
START_OBSTACLE_l 
END_OBSTACLE  1 
START_OBSTACLE_2 
END_OBSTACLE_2 


TIME_TYPE 
CCX)RD_TYPE 
COORD_TYPE 
COORD.TYPE 
COORD  TYPE 


08  08  1991  1000  00 

32U-PD004877 

32U-PD001870 

32U-PD004880 

32U_PD000885 


-  CALCULATE_MOE 

-  CALCULATE_MOE 

-  CALCULATE310E 

-  CALCULATE.MOE 

-  CALCULATE  MOE 


RUN_ID  wordt  gebmikt  om  de  naam  van  de  logfile  te  maken.  FIRST_ID  en  LAST_ID  bcpalen 
van  welke  simulaties  de  units  worden  weggeschreven  en  van  welke  simulaties  de  MOE’s  worden 
berekend.  START_TIME  is  van  belang  bij  bet  bcrekenen  van  de  DURATION.  De  start-  en 
eindcoordinaten  zijn  coordinaten  van  de  hindemissen  die  vcrdedigd  moeten  worden. 
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